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FORZWORD 


This  report  was  prepared  bj  the  Inglewood  Laboratory,  Aircraft  System 
Division  of  Thcapson  Products,  Inc,  under  DSAF  Contract  Ifo.  33(6l6)-3729, 
The  contract  was  initiated  under  Project  No.  3048,  "Aviation  Puela", 

Task  No,  30178,  "T?ffect  of  Fuel  Properties  on  Fuel  Systems",  It  was  ad¬ 
ministered  under  the  direction  of  the  Propulsion  Laboratory  and  l^terlals 
Central,  Directorate  of  Advanced  Systems  Technology,  Urlg^t  Air  Develop¬ 
ment  Division,  with  C.  R.  Hudson,  J,  H,  Johnson  and  Lt.  T.  L.  Oossage 
acting  as  project  engineers. 

This  report  oovers  the  work  performed  during  the  period  February  1955 
to  December  1958. 

L,  A,  Schlagel  was  the  project  manager  for  Thompson  Prodt2ct3,  Inc, 
and  g,  Flndl,  H,  Brande  and  H,  Edvards  conducted  the  experimental  work 
at  the  Inglewood  Laboratory,  Inglewood,  California, 
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ABSTRACT 


Ths  purpose  of  the  program  was  to  determine  certain  phjsical  and  therno- 
dynamle  propertiea  of  four  selected  aircraft  fuela.  These  fuels  were  MIL~P- 
5624c  Grade  JP-4,  a  MrL-F-25558  Grade  RJ-1»,  MII-F-25656  Grade  JP-6-H,  and 
Decalin,  a  commercial  hydrocarbon  mixture  of  the  els  and  trans  isomera  of 
decahydronapthalene . 

Specific  properties  evaluated  woret 

1.  Equilibrium  solubility  of  air,  nitrogen  and  ethane  in  the  JP-4 
and  RJ-l  fuels, 

2.  Effect  of  dissolved  air,  nitrogen,  and  ethane  on  the  viaeosity 
of  the  JP-4  and  RJ-l  fuels, 

3.  Effect  of  dissolved  ethane  on  density  of  the  JP-4  and  RJ-l 

fuels. 

4.  Effect  of  pressure  level,  agitation  rate,  and  rate  of  pressure 
change  on  the  evolution  rate  of  air  and  ethane  from  the  JP-4  and  RJ-l  ftials.  , 

5.  Therrflodynamlc  properties  of  each  of  the  four  teat  fuels  in¬ 
cluding  specific  heats,  enthalpy,  and  entropy. 

6.  Vaporization  characteristics  of  each  of  the  test  fuels  includ¬ 
ing  vapor  pressure,  equilibrium  vaporization  curves,  ainimua  reflux  curves, 
variation  of  vapor  molecular  weight  during  minimum  reflux  distillations,  and 
variation  of  density  during  minimum  reflux  distillations. 

7.  Variation  of  liquid  density  as  a  function  of  temperature. 

8.  Literature  survey  correlation  compilation  for  the  determination 
of  physical  and  thermodynamic  properties  of  petroleum  fhels. 

•Formerly  d63ignat«d  as  Shell  OMF,  Grade  C. 
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GKiOKi'JlL  li/TBCDUCriCN 


This  report  is  presented  in  fire  Parts: 
PART  I  -  KCtUILIERIUM  SOLUBILm. 


Tha  equilibrium  nolubillty  of  air,  nitrogen  and  ethane  and 
their  effect  on  iriacosity  were  doteminad  for  an  RJ-1  fuol. 

In  addition,  tho  equilibrium  solubility  of  ethane  and  ito 
effect  on  Tijcoaity  in  a  JP-4  fuel  were  determined.  Part  I 
Includes  gen8rall2ad  correlations  for  extension  of  the  data 
to  other  fuels  of  the  seme  typo.  Studies  of  the  equilibrium 
solubility  of  air,  nitrogen  and  carbon  dioxide  in  JP-4  are 
also  reported.  Thane  studies  performed  at  tho  Vest  Coast 
Laboratory  prior  to  this  program  were  reported  to  tha  Coordi¬ 
nating  Research  Council  at  an  earlier  date. 

PART  II  -  EVOLUTICN  iUTE  £?rUDlES. 

This  second  Part  ia  concerned  with  non -equilibrium  evolution 
rate  studies  on  JP-4  and  RJ-1  fuala.  The  studies  were  ac- 
ccEiplished  through  the  use  of  a  simulated  aircraft  fuel  tank. 
Rate  of  pressure  change,  agitation  rate  and  pressure  level 
are  the  variablea  evaluated.  It  became  increasingly  evident 
during  this  phase  of  the  program  that  the  phenanena  of  evo¬ 
lution  ia  extremely  cceplex.  While  tha  study  results  report¬ 
ed  are  informative,  it  has  not  yet  been  possible  to  correlate 
the  data  and  evolve  a  definite  statement  regarding  non- 
equilibrium  evolution.  The  tost  program  has  shown  however, 
that  certain  variables  are  probable  controllable  variables. 

PART  III  -  VAPCRIZATICN  CYCLE  AM)  THE  TiSRMCDYNAMIC  PR015RTIES  OP 
JP-4,  RJ-1  AJd)  JP-OrH. 

Certain  thermodynamic  and  physical  properties  of  each  of  ths 
test  fuels  are  described  in  Part  III,  Emphasis  was  placed  on 
the  conecruction  of  entropy-temperature  charts.  Tha  use  of 
these  charts  for  tho  study  and  evaluation  of  petroleum  frac¬ 
tions  of  the  JP  and  RJ  variety  is  discussed. 

PAi?r  IV  -  TAPCRIZAIICI  CYCL5  AM)  THSRMCDYNAMIC  HIOESRTIES  OF  DBCALIN. 


Certain  thermodynamic  and  physical  properties  of  Decalin  are 
described  in  Part  IV. 


Manuscript  reletissd  by  authors  September  i960  for  disbribution  ns  a  Technical 
Report, 
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PAfiT  V  -  GEftERALIZED  THERMODYNAMIC  AfJD  PHYSICAL  HIOPERTIES  OF  lETROIHUM 
TYPE  FUELS. 

Part  7  is  a  canpllatlon  of  rarious  phyaiochemlcal  correlatlona 
taken  frcm  aTallable  literature*  Conalderable  effort  went  into 
the  eraluation  of  a  great  Tariety  of  this  type  information  re¬ 
ported  in  technical  Journals*  The  charts  presented  in  Part  V 
are  believed  the  most  accurate  and  suitable  correlations  publish¬ 
ed  to  date  . 

This  Final  Technical  Report  is  treated  in  five  Parts*  A  general  Table  of 
Contents,  and  list  of  Illustrations  appears  at  the  beginning  of  the  report, 

A  list  of  Symbols  for  each  part  appears  at  the  end  of  that  part;  reference 
follow  each  part*  An  individual  set  of  calculations  end  an  Appeisdix  are 
also  included  where  applicable* 

The  overall  research  program*  as  initially  proposed*  was  designed  to  provide 
infoimation  for  various  aspects  of  the  SM64  Navaho  Program*  Numerous  change 
were  made  and  modifications  incorporated  during  the  life  of  the  program* 
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PAR.T  I  EQUILIBRIUM  SOLUBILITY 


INTRODUCTION 

An  engineering  report,  Ing.  Er.  No,  183,  was  published  by  Thompson 
Products  West  Coast  Laboratory  on  August  19,  1955.  This  report  was 
concerned  with  gas  solubilities  in  selected  military  fuels.  The 
studies  were  performed  under  an  Air  Force  contract.  After  the  pub¬ 
lication  of  this  report,  interest  in  a  new  fuel,  RJ-1’^,  and  a  pressur¬ 
izing  gas,  ethane,  arose.  An  investigation  of  gas  solubilities  in 
these  two  military  turbine  fuels  and  the  effect  on  viscosity  then  was 
conducted.  The  research  results  of  this  investigation,  conducted  under 
contract  No.  AF  33(616)3729,  are  presented  herein. 

OBJECTIVE 

The  investigation  was  made  for  the  following  purposes: 

1)  The  determination  of  equilibrium  solubility  of  ethane 

in  a  JP-4  fuel  at  temperatures  to  500*^  and  at  pressures 
to  400  paia, 

2}  The  determination  of  equilibrium  solubility  of  nitrogen, 
air,  and  ethane  in  RJ-1  fuel  at  temperatures  to  500“F 
and  at  pressures  to  400  psia,  and 

3)  The  determination  of  viscosities  of  both  JP-4  and  RJ-1 
fuels  as  functions  of  temperature  and  dissolved  gas 
concentrations. 

SUMM.\RY 

The  test  results  indicate  that  the  solubility  of  nitrogen  and  air,  in  the 
RJ-1  fuel,  increases  with  rising  temperature.  The  solubility  of  ethane, 
in  both  the  JP-4  and  RJ-1  fuels,  decreases  with  increasing  temperature. 
The  solubility  of  ethane  was  lower  in  the  RJ-1  fuels  than  in  the  JP-4 
fuels. 

*  Former  designation.  Shell  UMF,  Grade  C. 
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Apparent  irrejjularitjea  in  air  solubility  at  300“F  and  400°!F,  were 
obser/ed  in  both  /ucla.  These  Irregularities  are  believed  to  be  caused 
by  oxygen- consuming  chemical-reactions  such  as  gum  formation. 
Similar  irregularities  were  reported  previously.  See  Reference  1. 

METHOD  AND  PROCEDURE 


TEST  APPARATUS 

Figure  «  1  is  a  schematic  of  the  Teat  Apparatus  designed  and  built 

by  the  West  Coast  Laboratory  of  Thompson  Products,  Inc.  This  unit 
was  used  to  determine  liquid  viscosity,  gas  solubility  and  fuel  expansion. 
The  teat  cell  is  further  illustrated  in  Figure  2  . 

In  general,  the  methods  and  procedures  described  in  Reference  1  were 
followed.  Certain  small  changes  in  procedure  and  test  instrumentation 
were  incorporated  to  improve  accuracy. 

The  test  sample  was  contained  in  the  inner,  true-bore  glass  tube  of 
the  cell  as  shown  in  Figure  .  3.  This  tube  also  contained  a  glass 
ball  which  was  0.  0083  inch  smaller  in  diameter  than  the  tube.  Installed 
adjacent  to  the  tube  in  the  cell  wae  a  eteel  scale  which  indicated  the 
total  fuel  volume  at  test  conditions,  with  an  accuracy  of  iP.  3  per  cent. 

Encasing  the  inner  cell  was  an  outer  glass  tube.  The  heati.ng  fluid  was 
circulated  through  the  annular  space  bet-va  »n  the  two  tubes.  The  fluid 
was  a  high-flash-point,  water-white,  mineral  oil  to  which  a  flash 
inibitor  had  been  added. 

The  ends  of  both  tubes  of  the  test  cell  were  plugged  with  invar  metal 
end-caps  housing  "0"  ring  seals.  Two  "0"  rings  were  required  for  each 
end  of  the  inner  glass  tube.  A  suitable  compound  for  use  in  contact 
with  both  the  heating  oil  and  the  fuel  could  not  be  found.  Silicone 
Compound  No.  SE-360*  "O"  rings  were  used  in  contact  with  the  heat 
exchange  oil.  Silicone  Compound  No.  7180'^  ”0"  rings  were  used  at 
the  point  of  contact  with  the  fuel. 


■t-  Plastic  &  Rubber  Products,  Inc.  -  Identification  Numbers 
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The  end  caps  were  uupported  and  spaced  tjy  a  steel  sleeve.  The  sleeve 
was  mounted  on  pivots  so  that  the  cell  coulil  be  rocked  for  agitation  of 
the  fuel  and  rapid  dissolution  of  the  solute.  A  spirit  level  was  mounted 
on  the  sleeve  so  that  the  angle  of  the  tube  was  held  constant  within  i  30 
seconds  for  viscosity  measurement. 

Thermocouples  were  installed  in  each  end  of  the  test  assembly  to 
measure  temperatures  of  the  heating  oil,  the  liquid  fuel  sample,  and 
the  vapor.  Temperatures-. were  indicated  on  a  potentiometer  with  an 
accuracy  of  tl^F.  A  transducer  was  used  to  measure  the  pressure  in 
the  test  cell.  This  transducer  was'calibrated  within  psi.  To  reduce 
volume  external  to  the  test  cell,  the  transducer  cavity  was  filled  with 
mercury.  The  external  vjlume  v^as  less  than  0.  OZO. cubic  inch.  A 
schematic  diagram  of  the  test  installation  is  shown  in  Figure  1. 


Flexible  metal  hose  was  used  to  carry  heating  oil  to  the  test  cell  in 
order  to  allow  for  oscillation  of  the  assembly  during  test.  The  heating 
oil  reservoir  contained  thermostatically  controlled,  immersion-type 
electrical  heaters.  Also  installed  in  the  reservoir  was  a  water  coil 
which  was  used  for  cooling.  For  the  high  temperature  run,  additional 
electrical  bayonet-type  heaters  were  installed  in-line. 


The  gas  to  be  dissolved  was  contained  in  steel  bottles.  These 
reservoirs  w-ere  charged  with  pressurizing  gas  to  approximately  700 
psig.  They  were  immersed  in  an  oil  bath.  The  oil  bath  was  maintained 
at  a  constant  temperature  of  ISO^F  tl**.  A  Bourdon  tifce-type  pressure 
guage,  accurate  to  1 1  psi,  and  a  differential  pressure  transducer 
calibrated  to  ilpsi  were  used  to  measure  the  pressure  in  the  steel 
bottles.  One  reservoir  was  used  as  a  reference  in  order  to  utilize  a 
d iff e r e nt ia  1  —  p r e s s u r e  transducer  with  a  0  to  ZOO  psi  range.  To  increase 
the  accuracy  of  tne  test  results,  this  transducer  was  preferred  to  a 
conventional  transducer  or  pressure  guage  having  a  0  to  700  psi  range. 

Connections  to  the  test  cell  verc  made  by  flexible  capillary  tubing. 


The  cell,  shown  in  Figure  4,.  was  used  to  measure  the  gas-free 
vapor-pressure  of  the  fuel.  The  reflux  column  was  cooled  by  a  dry  ice 
and  acetone  bath.  The  cell  could  be  sealed  from  the  reflux  column.  A 
hot  oil  bath  was  provided  for  temperature  control  of  the  cell  assembly 
as  shown.  Cell  pressure  was  measured  with  the  calibrated  transducer 
described  above.  Temperatures  were  measured  by  means  of  a  thermo¬ 
couple  (and  potentiometer). 
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TEST  PROCEDURE 

Each  teat  fuel  waa  subjected  to  the  following  physical  and  chemical  in- 
apectiona  in  accordance  with  the  applicable  milila.  y  opccification. 

Diatillation  --  per  VV- L-79!,  Method  1001,7 

Gravity  -  per  VV-L-79I,  Method  401.  3 

Reid  Vapor  Preaeure  -  per  VV-L-79!,  Method  1201,  4 

Initial  oxygen  and  nitrogen  contents  of  the  fuels  were  determined  by  a 
special  laboratory  apparatus  which  is  shown  in  Figure  5.  The  fuels 
were  degassed  by  means  of  mercury  pumps  and  the  gases  were  stored 
in  separate  burettes.  After  complete  evolution  of  the  gases,  Orsat 
analyses  were  made  on  the  total  gas  sample  collected. 

Solutes  of  high  purity  were  used  in  the  investigation.  Analyses  of  each 
solute  are  presented  in  Table  I. 

The  volume  of  the  test  cell  and  the  gas  reservoir  were  determined 
independently  by  reduction  of  presstire  in  each  cell  to  one-half  atmos¬ 
phere.  The  volume  of  gas  removed  was  measured.  Inasmuch  as  the 
cells  and  the  gas  reservoirs  were  maint<u  ned  at  a  constant  temperature 
the  volume  of  gas  removed  was  a  direct  measurement  of  the  cell  volume. 

The  volume  increase  of  the  tost  cell  due  to  thermal  expansion  at  500°F 
was  found  to  be  less  than  0.  3  per  cent.  This  change  in  volume  did  not 
significantly  affect  the  calculations. 

The  viscosity  calibrations  were  reported  provioi.sly.  See  Reference  1. 

Solubility  and  viscosity  measurtin  enis  were  made  in  thf'  test  cell  hv 
use  of  the  following  procedure: 

Before  each  run  the  test  cell  and  b^^ll  were  cleaned  witn  a 
sulfuric  acid-potassium  dichromate,  glass-cleaning  solution. 
When  required,  the  end  caps  were  cleaned  with  hydrochloric 
acid.  All  parts  were  thoroughly  rinsed  with  water  and  acetone 
and  dried  before  assembly. 
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The  outer  glass  tube,  and  the  inner  glass  tube  containing  the 
glass  ball,  were  assembled  inside  the  metal  spacer  with  one 
of  two  end  caps.  At  a  known  temperature  a  meastired  volume 
of  the  fuel  to  be  tested  was  placed  in  the  cell.  The  quantity  was 
adjusted  so  that  the  vapor-liquid  ratio  at  the  run  temperature 
would  be  approximately  0.  15.  The  second  end  cap  then  was 
fastened  in  place.  The  cell  was  placed  on  a  test  stand  under  a 
vapor  hood.  The  capillary  gas  line,  the  heating  fluid  lines, 
and  the  thermocouple  leads  were  then  connected.  The  fuel 
temperature  and  lic|uid  level  in  the  cell  were  recorded.  The 
liquid  level  was  indicated  by  a  metal  scale  inside  the  cell,  read 
with  the  cell  in  a  vertical  position. 

The  heating-oil  circulating-pump  and  the  heaters  were  then 
turned  on  and  the  cell  oscillated  to  agitate  the  fuel.  Run 
temperatures  were  reached  in  10  to  50  minutes. 

As  soo.n  as  the  cell  was  stabilized  at  the  run  temperature 
another  liquid  level  reading  was  taken  and  the  cell  press..re  and 
temperature  were  recorded.  The  cell  was  then  positioned  at  an 
angle  of  approximately  15°  above  horizontal  for  the  viscosity 
determination;  this  angle  was  constant  within  ±30  seconds  for 
all  measurements. 

After  this  measurement  the  cell  was  again  set  into  oscillation. 

A  quantity  of  gas  was  metered  through  the  capillary-line  valve 
into  the  cell.  The  flow'-rate  was  adjusted  so  that  the  cell  pressure 
never  exceeded  the  chosen  pressure  point  by  more  than  5  psi. 
Agitation  was  continued  for  10  mir.atec  after  the  cell  pressure  had 
stabilized.  The  liquid  level,  roll  time,  cell  pressure,  and  fbe 
pressure  drop  in  the  gas  reservoir  were  recorded. 

Additional  gas  was  introduced  into  the  cell  and  the  above  process 
repeated  until  data,  at  a  minimum  of  for.r  pressure-levels,  were 
recorded.  T!ie  highest  pressure-level  investigated  was  415  psia. 
This  procedure  was  repeated  at  various  fuel  temperatures  to 
500°F. 

The  data  obtained  from  the  above  procedure  wore  reduced  in  the 
followdng  man.ncr. 

1.  Roll  time  was  converted  to  absolute  viscosity  by  use  of  a 
caiibration  curve. 


2.  The  calculatio.ns  used  in  the  reciuc:tion  of  the  sol'..bility  data 
are  described  in  the  3oli.''ilitv  ralrulatioas. 

".ADD  TH  6C-7t7  7 


SOLUBILITY  CALCULATIONS 


The  solubility  calculations  were  based  on  Dalton's  law  of  partial 
pressures;  Each  constituent  contributes  to  the  total  pressure  a  partial 
pressure  equal  to  the  pressure  it  would  exert  if  it  alone  were  present 
aj  the  given  temperature  in  the  volume  occupied  by  the  mixture. 
Reference  3.  From  the  meager  quantity  of  mixture  data  available  in 
the  literature,  Dalton's  law  appears  to  be  quite  satisfactory  for  mixtures 
at, pressures  below  50  atmospheres. 

Dalton 's  law  can  be  expressed  mathematically  as  follows; 


=  £pi  =  Pi 


+  P2  +  P3 


Pi  ''■i  =  Pi  yi 


Pi  = 


(  Pi  Vj 


By  definition  Z  =  or  similarly  Z-  -  ^ _ L 

RT  ‘  RT 


and  Z. 


P  V 
t  t 


'i  I  Zj 
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A  simplified  block  diagram,  as  shov/a  below,  will  be  used  to  illustrate 
the  relationships  derived  from  Dalton's  law.  These  relationships  were 
used  to  determine  the  solubilities  of  the  three  solutes  {air,  nitrogen 
and  ethane  )  in  the  two  solvents,  RJ-1  and  JP-4. 


It  is  evident  that  once  Valves  2  and  3  have  been  closed,  all  gas  leaving 
Reservoir  A  must  enter  the  test  cell.  Reservoir  B  was  used  only  to 
provide  a  constant  reference  pressure. 
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The  quantity  of  gas  entering  the  cell  can  be  expressed  in  terms  of  the 
compressibility  forrh  of  the  equation  of  state, 

■  AR  V. 

An  -  *  *  _  An- 

Z  RT,  “  (5) 

r  r 


In  order  to  determine  the  quantity  of  pressurizing  gas  that  has  been 
dissolved,  all  that  is  necessary  is  to  determine  the  quantity  of  gas  that 
remains  in  the  ullage.  The  difference  between  the  quantity  of  pressur¬ 
izing  gas  which  enters  the  cell  and  that  which  remains  in  the  ullage, 
is  the  quantity  added  to  the  solution- 

Ud  =  A^r  -  "v,x 

To  evaluate  n^  Dalton's  law  is  utilized, 
vx. 

The  equation 


(7) 


which  states  that  the  total  cell  pressure  is  due  to  the  combined  partial 
pressures  of  the  pressurizing  gas,  the  residual-air  partial  pressure, 
and  the  fuel-vapor  pressure.  The  fuel-vapor  pres3\;re  and  air-partial 
pressure  w'as  assumed  to  be  a  function  of  temperature  alone, 

Px  "  ■  'Pa  *  Pf  >  (8) 


ThejiCalculation  of  pj^  can  be  made  when  the  cell  pressure  prior  to  the 
addition  of  the  pressurizing  gas  is  known. 
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The  quantity  of  preaeurizing  gas  in  the  ullage  can  be  determined 
utilizing 


PxVv 

“v  ^  Zx  RTv 


combining  equations  (5),  (6)  and  (9}> 


Pj  Vi 

Zv  RTi 


(9) 


(10) 


The  Ostwald  coefficient  is  used  to  put  Equation  (10)  into  a  more  useful 
form.  The  Ostwald  coefficient  is  defined  as  the  volume  of  solute, 
measured  at  the  temperature  and  pressure  at  which  it  is  dissolved, 
per  unit  volume  of  solvent. 


Hd  Zx  RTy 
Px 


(12) 


Ox  = 


M  Zx  RTy 
PxVl 


(13) 


Another  useful  form  to  express  the  quantity  of  Equation  (1)  is  to  state 
the  volume  dissolved,  per  unit  weight  of  solvent,  measured  at  standard 
conditions. 
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n  ,  Z  RT‘ 
d  X 


=  Z’  V, 


n_,  Z  RT' 

d  X 


"When  air  and  nitrogen  are  used  as  pressurizing  agents,  certain  minor 
modifications  of  the  foregoing  equations  are  required.  Dissolved  air 
has  been  arbitrarily  defined  as  the  combination  of  nitrogen  and  oxygen 
dissolved  in  the  solvent,  not  necessarily  in  the  same  ratio  as  found  in 
the  atmosphere.  This  definition  is  necessary  due  to  the  somewhat 
greater  solubility  of  oxygen,  as  compared  to  nitrogen,  in  hydrocarbon 
fuels. 

When  air  is  used  as  the  pressurizing  agent.  Equation  (8)  simplifies  to 


P  P 

t  c  -  F 


Since  air  is  already  dissolved  in  the  e  v,vent.  Equation  (6)  must  be 
modified  to  include  the  gas  initially  dissolved.  The  total  gas  in  the 
cell  at  any  time  is  defined  by 


n  +  n  , 
V  d 


At  initial  conditions  prior  to  pressurization  from  the  reservoir  this 
becomes 

*^t  initial  *^y  ^  ''i 
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After  gaa  haa  been  added  from  the  reservoir  the  total  air  in  the  cell 
is  the  eum  of  that  initially  there  plus  that  added  from  the  reservoir: 


t  initial 


+An  -  n  +  n,  + 
r  y  1 


(20) 


Combining  Equation  (13)  and  Equation  (20) 


n 


n. 


"v  "  "y  ”i  “v 


When  nitrogen  is  used  as  the  pressurizing  agent  a  different  correction 
to  Equation  (8)  is  required.  This  is  due  to  the  fact  that  the  air 
initially  in  the  cell  ullage  consists  of  approximately  21%  oxygen  and 
79%  nitrogen.  Thus  for  nitrogen  calculations.  Equation  (8)  becomes 


n 


^c 


(Pp  +  0.21  p^) 


(22) 


To  summarize,  the  relationships  used  for  air,  nitrogen  and  ethane 
solubilities  are  as  follows: 

Air 


P  -  ^X/Vx  V  n  V® 

rl  r2\(  r  \  ^  14.7  v,  initial^  14,7  1,  a 


Zrl  Zr2/\T, 


^v  initial 


(P  ,  P„)  V 

'  t,  c  -  F  V 

- 27 - tT: - 


Z  T 

a  V 


(^^t.c  -PF)Vl 


(23) 
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v«  “.T"  [Ptc-<>’F  t  Pa»] 

W  =  Tv  ^p“  ~  (28) 

The  compressibility  factors  of  air  and  nitrogen  are  nearly  unity  in  the 
range  investigated  and  were  neglected  to  simplify  the  calculation. 

Ethane  compressibility  factors  were  computed  from  d^ta  in  References 
2,  3,  4,  and  5.  The  information  contained  in  these  References  was 
cross-plotted  and  converted  to  a  form  of  the  virial  equation  of  state, 
i.  e.  ,Z«lf  B(T)  P+  C(T)P^  +  D{T)P^.  Figure  6  shows  the  results 
of  the  virial  equation,  plotted  at  the  test  temperatures.  Figures  7 
through  9  show  the  values  of  the  second,  third  and  fourth  virial 
coefficients  used. 


•ADD  XR  60-767 


15 


DISCUSSION 

The  test  results  demonstrate  the  general  trend  indicated  in  Reference  I, 
i.  e. ,  solubility  decreases  with  increasing  molecular  weight  of  solvent 
fuel. 

As  anticipated,  ethane  is  extremely  soluble  in  the  test  fuels.  Solubility 
of  this  solute,  as  determined  by  the  Ostwald  Coefficient,  is  pressure 
dependent.  This  is  due  to  the  large  deviation  of  the  ethane  compressi¬ 
bility  factor  from  unity  and  the  radical  change  in  liquid  density  caused 
by  this  solute.  The  deviation  can  best  be  expressed  mathematically. 

■  '■x  ’'l  '  '*'1,  \  •’x 


Differentiation  of  Equation  (1)  with  respect  to  p  at  constant 
temperature  results  in 


O  = 

X 
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Ideally, 


!  3 


S>  z 


ap 


=  a  constant. 


BP 


ap 


s  p. 


A 


0  and  thus 


the  change  of  Ostwald  Coefficient  with  pressure  is  zero.  In  the  case 
of  the  ethane-Fuel  A  (JP-4)  and  the  ethane  -  Fuel  B  (RJ-1)  systems. 


^  and  /  3/^l' 


do  not  equal  zero.  Therefore,  the 


^  Px  Vt 


Ostwald  Coefficient  can  increase  or  decrease  with  pressure.  The  slope, 
[  ^  ^x  ,  is  dependent  upon  the  rejat^nsh'p  \etween  »  Z^ 


1.  e. 


and  p  .  Both  Fuel  A  and  Fuel  D  show  deviations  frorri  ideality.  See 


Figures  12  and  15.  In  the  low  temperature  region  (100“F  to  2UU'^F), 


!  ^  ^x  ^  is  negative.  At  the  high  temperatures,  (300‘T*  to  500“F) 

V  ^Px  /  T 


X 

a  o 


is  positive,  indicating  that  in  the  200°F  to  300'T'  range 


ap. 


/  ao 


I- 


3p. 


=  0:  That  is,  a  point  is  re^tChed  where  the  system  behav’es 


X  /  T 
ideally. 


.Air  and  nitrogen  behave  ideally  in  both  JP-4  and  RJ-1  fuels,  which 

/  5^1,  \ 

could  be  predicted  since!—;;- -  i  and  I - i  are  both  t-t-  triv 

V  ®  '’x  /t  \  3  ’’x  jr 


eijual  to  zero.  In  the  200"  F  to  400"F  range  the  apparent  solub.llt  f 
air  increases  somewhat  rapidly.  This  is  due  to  a  visible  c  hemic.n 
reaction  (presumably  oxidation)  which  causes  formation  of  de[iosits.  A 
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true  eolution  effect  is  not  achieved.  The  400^^  runs  on  nitrogen  do  not 

AOx 

follow  the  — —  trend  shown  by  the  runs  performed  at  other  tempera¬ 
tures,  An  explanation  for  this  behavior  has  not  been  deduced. 

The  solubility  of  various  pressurizing  gase's  present  in  the  test  fuels 
were  compared  and  are  illustrated  in  Figures  10  and  11.  Specific 
information  on  carbon  dioxide,  air  and  nitrogen  in  JP-4  fuels  was 
obtained  from  Reference  1.  The  data  has  been  extrapolated  to  the 
critical  temperature,  since  by  definition,  the  Ostwald  Coefficient  must 
equal  1.  0  at  the  critical  temperature  of  the  solvent. 

There  does  not  seem  to  be  any  theoretical  reason  for  the  apparent 
reversal  of  solubility  of  ethane,  in  Fuel  B,  with  increasing  temperature. 
See  Figure  11.  Due  to  the  method  used  in  the  evaluation  of 
solubilities,  the  ethane  data  at  the  higher  temper  atures  is  believed  to 
be  accurate  only  to  ilOfa.  ♦  This  deviation  is  a  result  of  the  compara¬ 
tively  large  changes  in  compressibility  factor  resulting  from  changes  in 
temperature  and  pressure.  With  this  deviation  ,  it  is  possible  that  the 
Ostwald  Coefficient  does  not  fall  below  1.  0  and  a  solubility  reversal  does 
not  occur. 


*  The  compressibility  data  is  believed  to  be  accurate  to  0.  5%, 

temperature  instrumentation  to  ±1.  0%  and  pressure  instrumenta¬ 
tion  accurate  to  ±1.  0%.  However,  a  1%  temperature  deviation 
could  mean  a  change  of  as  much  as  1.  S'F  in  the  reservoirs.  This 
change  could  increase  the  compressibility  factor  deviation.  In 
the  calculations,  any  error  in  compressibility  factor  is  multiplied 
due  to  the  mathematical  methods  employed.  Thus,  from  relative¬ 
ly  small  instrumentation  errors,  large  errors  in  results  are 
poseible. 
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Air  and  nitrogen  have  compreesibilities  nearly  equal  to  Unity  and  are 
relatively  inaensitive  to  comprcsflibility  difficulties.  The  data  for 
these  gases  should  approach  an  accuracy  of  ±2%. 

The  solubility  data  obtained  in  terms  of  the  Ostwald  Coefficient  are 
shown  in  Figures  12  through  15,  Figures  16  through  19 
show  solubility,  in  terms  of  weight  of  solute  dissolved  per  weight  of 
solvent,  as  a  function  of  pressure. 

The  effect  of  dissolved  gases  on  the  viscosity  of  the  test  solvents  Is 
shown  in  Figures  20  through  23.  The  ability  of  ethane  to 
lower  viscosity  is  much  more  pronounced  at  lower-liquid  temperatures. 

Figures  24  and  25  show  the  effect  of  dissolved  ethane  on  the 
density  of  the  two  test  fuels,  JP-4  and  RJ-1, 
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SYMBOLS 


D, 

m 

Constant* 

n 

m 

Moles  of  pressurising  gas 

0 

m 

Ostwald  coetficient 

P 

m 

Pressure 

P 

w 

Partial  pressure 

R 

m 

Gas  constant 

T 

M 

Temperature 

V 

m 

Volume 

W 

m 

Weight 

Z 

m 

Compressibility  factor 

A 

m 

Change 

m 

Solvent  density  at  charge 
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SUBSCRIPTS  AND  SUPERSCRIPTS 
Air 

c  =  Cell 

d  =  Solute  diasolved 

e  s  Ethane 

£  =  Fuel 

i  =  Solute  initially  dissolved  in  solvent 

L  =  Liquid  solvent 

N  =  Nitrogen 

r  =  Reservoir 

t  =  Total 

V  =  Cell  ullage  or  vapor  space 

X  =  Pressurizing  gas 

y  =  Solute  initially  in  vapor  space 

z  =  Last  constituent 

o  =  Standard  conditions 

1  , 2,  3  =  Constituents  1,  2,  3,  etc. 
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Figurn  No.  2  «  Photograph  of  Teat  Cell  Aafembly 


-7 


No.  1!  -  Ostwald  Solubility  Coofficifiit  of  V.irioiia  Gases  in  R J- 1  Fuel. 
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Nv'.  -  (.^hIw.iIcI  Solubility  Coi'ffiriont  of  Ethane  in  RJ-1  Fuel. 


an« 


afINiird 


Solubility  of  Nitrogen 


Solubility  of  Ethane  in  Fuel  B  -  RJ-l. 
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ect  of  Diasolved  Ethane  on  Density,  Fuel  j 


No.  25  -  Effect  of  Di**olved  Ethane  on  Deneity,  Fuel  B 


_ SOLUTE  ANALYSES 

1.  AIR 

H2O 
CO2 

Illuminanta 

O2 

CO 

Hydrocarbons 

Nz 

2.  NITROGEN 

Purity  {minimum) 

3.  ETHANE 


Research  Grade 


PART  II  EVOLUTION  RATE  STUDIES 

INTRODUCTION 


This  Part  ia  one  of  a  series  of  interioi  reportst  each  of  which 
is  offered  in  partial  fulfillment  of  Research  Program  AF  33(616)3729. 
This  program  was  initiated  in  June,  1956.  The  results  reported  herein 
cover  Phase  4  of  the  Research  Program.  Phase  4  represents  a  study  of 
the  evolution  rates  of  three  pressurizing  gases  from  two  selected, 
military  fuels.  The  two  fuels  are;  Fuel  A,  (JP-4)  and  Fuel  B  (RJ-1). 


OBJECTIVE 

The  purpose  of  Phase  4  of  the  Research  Program  was  to  determine  the 
effect  of  fuel  type^plus  pressure-change  and  agitation  ratesj^on  the 
evolution  rate  of  three  pressurizing  gases.  Test  conditions  simulated, 
as  closely  as  practical,  conditions  occurring  in  present  aircraft  and 
missiles.  Conditions  anticipated  in  future  aircraft  and  missiles  were 
also  simulated  and  tests  performed  in  this  environment. 

SUMMAR  Y 


Experimental  data  were  obtained  for  the  evolution  rates  of  four 
systems; 

Air  -  JP-4  Air  -  RJ-1 

Ethane  -  JP-4  Nitrogen  -  RJ-1 

These  tests  determined  that  the  evolution  rate  of  ethane  is  consider¬ 
ably  higl-er  than  that  of  air.  Further,  agitation  level  has  a  pronounced 
effect  on  evolution  rate.  The  effect  of  pressure  level  on  evolution 
rate  was  not  observable  in  the  air-fuel  systems,  Imt  was  conspicuous 
in  the  ethane -fuel  system. 
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Comparison  between  the  test  data  and  the  work  of  Schweitzer  and 
Szebehely  *  showed  little  correlation.  This  is  believed  due,  at  least 
in  part,  to  different  experimental  approaches  and  control  of  variables. 


TEST  EQUIPMENT  AND  PROCEDURE 


The  equipment  of  this  phase  of  the  program  was  designed  to  simulate 
an  aircraft  fuel  tank.  Certain  environmental  test  conditions  were  then 
imposed  upon  this  tank.  The  conditions  approximated  those  occurring 
in  high  performance  missiles  and  aircraft. 

Initially,  it  was  intended  to  schedule  a  programmed  flight  plan  through 
take-off  and  cruise,  monitoring  both  evolution  rate  and  tank,  ullage 
composition.  One  of  the  environmental  conditions  to  he  imposed  on  the 
tank  during  the  cruise  portion  of  this  simulated  flight  was  heat.  Skin 
temperatures  to  1000*F  and  pressures  to  200  psia  were  to  be  investigated. 
The  tank  ullage  composition  was  to  be  monitored  using  a  rapid- scan  mass 
spectrometer.  Due  to  time  and  monetary  considerations,  this  portion  of 
the  program  was  deleted. 

The  test  tank,  however,  was  constructed  to  meet  the  initial,  severe 
environmental  conditions.  The  tank  is  contained  in  a  stainless- steel 
oven  where  radiant  heaters  (stainless-steel-sheath)  supply  the  thermal 
energy.  A  special  steel  teat  chamber  was  built  to  house  the  entire  unit. 
See  Figures  26  through  2?  for  additional  detail. 

A  considerable  effort  went  into  the  development  of  a  suitable  means  of 
providing  agitation.  There  appears  to  be  no  set  standard  of  agitation, 
nor  is  there  any  approved  test  method  for  simulating  aircraft  agitation. 
Initially  it  was  proposed  to  develop  a  stirrer-type  agitator  that  did 
not  employ  seals.  Tank  leakage  could  not  be  tolerated  because  of  the 


*  Schweitzer,  P.  H.  and  Szebehely,  V.  C.  Gas  Evolution  in  Liquids 
and  Cavitation.  Journal  of  Applied  Physics.  Volume  21.  December, 
1950.  pp  1218  -  1224. 
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radiant  heater*,  even  though  provision*  for  inerting  the  oven  space  were 
provided. 

The  agitator  that  wa*  developed  wa*  a  propeller  type.  See  Figure  ^0, 

The  propeller  wa*  driven  by  an  eccentric  shaft  and  pinion  device  en¬ 
closed  in  a  bellows.  The  entire  construction  wa*  of  stainless  steel. 
Trouble  developed  however,  when  the  bellows  failed  repeatedly  after  a 
short  period  of  operation. 

A  second  type  of  agitator  was  then  developed.  This  agitator  consisted 
of  a  rocker  mechanism  which  tilted  the  oven  and  tank  about  its  center 
of  gravity.  See  Schematic  shown  in  Figure  2S*  The  amplitude  of 
deflection,  measured  at  the  end  of  the  oven,  is  ^  inch.  The  agitation 
rate  is  varied  by  changing  rocking  frequency, by  using  different  pulley 
ratios,  while  using  a  constant  speed,  one-half  HP  motor.  Four  agitation 
frequencies  were  arbitrarily  chosen:  350  CPM,  240  CPM,  150  CPM  and 
zero  CPM.  This  method  q£  agitation  proved  feasible  and  was  used  to  • 
obtain  the  data  presented  in  this  Supplement. 

A  35  liter  aircraft-oxygen-bottle  was  used  as  the  vapor  reservoir.  The 
bottle  wa*  connected  to  the  tank  with  a  flexible  tubing  line  to  damp  the 
tank  vibrations.  A  sump  container  was  inserted  between  the  reservoir 
and  the  tank  to  prevent  any  liquid  entrainment  losses  from  entering  die 
reservoir. 

Tank  and  reservoir  pressures  were  recorded  using  pressure  transducers 
and  a  light  beam  oscillograph.  Tank  and  reservoir  pressures  were 
monitored  visually  by  means  of  pressure  guages.  Temperatures  were 
measured  usingiron-constanfen  thermocouples  and  galvanometer  read¬ 
out  instruments.  All  temperature  measurements  were  recorded  manually. 

The  test-procedure  steps  follow; 

1.  Flush  tank  with  fuel  to  remove  possible  contamination. 

2.  Fill  tank  with  fuel  to  prescribed  level. 
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3.  Evacuate  reservoir. 

4.  Record  initial  conditions. 

5.  Pressurize  tank  to  desired  level.  Agitate  with  rocker  to 
speed  dissolution  of  pressurizing  gas. 

6.  Measure  and  record  the  tank  liquid  and  vapor  temperatures, 
plus  the  reservoir  temperature,  after  equilibrium  is 
reached  in  tank. 

7.  Set  tank  agitation  rate  by  pulley  adjustment. 

8.  Turn  on  recording  instrumentation. 

9.  Vent  into  reservoir  to  drop  tank  pressure.  Rate  of  pressure 
drop  is  manually  controlled  for  a  linear  decline. 

10.  Close  vent  valve  when  desired  final  pressure  is  reached. 

11.  Turn  off  recording  instrumentation  and  note  temperatures. 

12.  Repeat  steps  6  through  11,  over  the  desired  pressure  range, 
after  equilibrium  is  reached  in  the  tank. 

Static  equilibrium  was  evaluated  by  noting  pressure  changes.  Initially, 
if  a  pressure  change  was  not  noted  within  a  lO-minute  period,  tank 
equilibrium  was  assumed.  This  time  period  proved  to  be  insufficient 
and  the  procedure  was  modified  to  allow  a  minimum  of  |-hour  for  the 
system  to  come  to  equilibrium  after  a  test  run.  Where  air  and  nitrogen 
vere  the  pressurizing  gases,  a  ^-hour  period  was  satisfactory.  Ethane 
required  much  longer  periods  of  time,  ranging  from  one  to  six  hours. 
This  time  period  was  dependent  to  a  large  degree  on  the  previous  rale  of 
pressure  change.  If  the  pressure  change  rate  had  been  rapid,  a  longer 
period  of  time  would  have  been  required  for  the  ethane  fuel  system  to 
achieve  equilibrium. 
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E  VOLUTION  RATE  CALCULATIONS 

The  evolution  rate  calculcTtions  were  based  on  a  system  weight  balance. 
The  system  consisted  of  a  fuel  tank  and  an  evacuated  chamber  called 
the  reservoir.  See  sketch  below. 


Note,  from  the  above  sketch,  that  when  the  vacuum  system  valve  and 
pressurizing  system  valves  have  been  closed  any  gas  or  vapor  that 
leaves  the  tank  must  enter  the  reservoir.  Expressed  mathematically 


Anr 


AP  V 
“r  *  r 

Z,  RT 
r  r 


AP^V^ 
Zy  RT“ 


Ane 


For  the  purpose  of  these  calculations,  Z^  is  assumed  to  equal  1.  0  at 

pressures  below  atmospheric.  This  is  approximately  true  for  both  the 

fuel  vapors  and  the  pressurizing  gases  used  in  this  study.  For  this 

reason  Z  has  been  omitted  in  all  calculations, 
r 


MALD  TH  60-767 


53 


AP 

The  uie  of  a  linear  change  rate,  i.  e,  *  K  gives  the  evolution  rate 

for  the  total  fuel  charged.  . 


The  weight  of  fuel  charged  can  be  determined,  when  the  volume  of  fuel 
charged  and  its  density  at  the  charge  temperature  are  known. 


W  » 


The  evolution  rate,  to  be  of  moat  use,  should  be  based  on  a  unit  quantity 
of  fuel.  Therefore  the  evolution  rate  per  unit  weight  of  fuel  can  be  ex¬ 
pressed  as  follows: 


Due  to  practical  difficulties,  such  as  changes  in  density  of  the  liquid 
fuel  due  to  temperature  and  pressurizing  gas,  certain  corrections  are 
required.  In  the  cases  of  air  and  nitrogen,  the  effect  of  the  pressuriz¬ 
ing  gas  on  liquid  density  is  negligible.  In  the  case  of  ethane  the  effect 
is  greater  and  the  correction  required  is  of  a,ppreciable  magnitude. 

The  fuel  tank  ullage  was  kept  small,  i.e. ,  small  vapor-liquid-ratio, 
during  all  tests  to  maximize  the  effect  of  gas  evolution.  Therefore  small 
temperature  changes,  although  causing  small  changes  in  fuel  volume, 
caused  relatively  large  changes  in  ullage  volume.  Using  Bureau  of 
Standards  reference  data*  for  the  volumetric  expansion  coefficient,  a 


♦  Cragoe,  C.  S.  Bureau  of  Standards  Misc.  Publication  No.  97.  1929- 
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calculation  of  the  change  in  liquld-fuel-volume  for  a  temperature 
change  can  be  made.  This  calculation  provides  a  means  for  computing 
the  change  in  ullage  volume. 


Val  =  V,  OL  (Tl  -  T,) 


Vl  =  ["n-OCtTL-Ti)] 

n  '  V.  .  .  Vj  -  Vi  r  ^  cC(Xl-  T;)! 

Substituting  this  correction  back  into  Equation  (8) 


£  = 


pPr  Vr 

L 


APy  (  Vs  -  Vj  [T»  CSCiTh-  T^l 


Zv  Ty 


Equation  (12)  was  used  for  computing  the  evolution  rates  of  air*  in  Fuel 
B  (RJ-l).  An  additional  correction  for  vapor  pressure  is  required  in 
the  case  of  Fuel  A  (JP-4).  The  correct. on  factor  for  Fuel  A  is  de¬ 
veloped  later. 


*  The  term  air  as  used  here  denotes  the  mixture  of  oxygen  and  nitrogen 
leaving  the  tank  and  entering  the  reservoir.  Due  to  the  somewhat 
selective  solubility  of  the  fuels,  i.  e.  ,  oxygen  being  more  soluble 
than  nitrogen,  the  actual  gas  evolved  from  the  liquid  will  be  somewhat 
oxygen  enriched  when  compared  to  the  gas  in  equilibrium  in  the 
vapor  space.  The  deviations  are  small  however,  and  for  most 
purposes  negligible. 


*  Cragoe,  C.  S.  Bureau  of  Standards  Mlac.  Publication  No*  97* 
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Due  to  the  very  high  solubility  of  ethane  in  petroleum  fuels,  the  liquid 
fuel  densities  deviate  considerably  with  pressurizing- gas  pressure- 
level.  The  density-temperature  effect  is  also  amplified.  This  lajtcr 
problem  is  due  to  the  fact  that  the  solubility  of  ethane,  in  petroleum 
fuels,  increases  quite  rapidly  when  temperatures  decrease  from 
ambient.  Figures  3I  and  32  show  the  density  variation,  due  to 
pressurizing- gas  pressure  at  various  temperature  levels,  when  equilib¬ 
rium  conditions  are  maintained, 

A  calculation  of  the  liquid  and  ullage  volumes  at  any  pressure  level  can 
be  made  using  this  data  and  interpolating  as  required.  This  is  done 
using  the  weight  of  fuel  charged,  from  Equation  (7),  which  remains 

constant;  and,  the 
condition. 


w  =  /^i  Vi  _  (13) 


W 

sic  definition  of  density,  i.  e. 


Since  the value  is  at  run-temperature  condition,  further  temperature 
correction  (i.  e.  ,  expansion  coefficient  correction)  is  not  required. 

In  actual  practice  equilibrium  conditions  are  not  obtained  diring  a  test 
run  except  at  the  initial  starling  point,  and  prior  to  a  decrease  in 
pressure.  The  deviation  from  equilibrium  is  a  function  of  many  vari¬ 
ables.  A  few  of  these  variables  are;  rate  of  pressure  change,  surface 
area,  viscosity  and  nucleation  sites.  To  accurately  obtain  the  dynamic 
densities  required  for  a  rigorous  calculation,  quite  elaborate  equipment 
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FUEL  DENSITY 


would  be  required.  As  correction  for  deviation  from  equilibrium  is 
just  one  of  the  correction  factors,  useful  only  on  the  ethane  pressuriz 
ing  system,  it  was  felt  that  approximations  were  in  order.  The 
approximation  consists  of  using  the  arithmetic  mean  density  from 
Figure  31  and/or  Figure  32  when  both  the  initial  and  final 
equilibrium  pressures  are  known.  (Where  the  final  pressure  is  that 
pressure  which  the  system  finally  attains  after  a  pressure  change). 
The  arithmetic  mean  was  used  for  its  simplicity  and  the  fact  that  in 
the  pressure  increments  used,  the  density  deviations  are  nearly  a 
linear  function. 


r  p 

'  1 _ 

Pui^  Pu 

v.  ^ 


Obviously  under  conditions  of  a  pressure  change  at  equilibrium,  this 
approximation  becomes  nearly  exact.  Under  conditions  of  non¬ 
equilibrium  the  correction  error  increases  to  a  maximum.  This 
maximum  error  should  not  exceed  more  than  an  estimated  10  per  cent 
of  the  change  in  density  nor  less  than  1  per  cent  of  the  average  density 
The  diagram  below  may  clarify  this  statement. 


ETH.^NE  P.ARTLAL  PRESSURE 


■0-ru, 


The  curve  E-H  represents  a  typical  density-pressure  plot  similar  to 
Figures  and  32  obtained  under  equilibrium  conditions.  Curve 

D-H  and  D'-H  represent  the  plots  obtained  by  venting  the  tank  from 
pressure  level  "C"  to  pressure  level  "A"  under  non-equilibrium  con¬ 
ditions  (D'-H  being  more  rapid  than  D-H).  Curve  D-F  represents  the 
change  in  density  after  the  vent  is  shut  and  equilibrium  conditions  arc 
obtained.  The  average  density  during  a  test  run  would  be  that  deter¬ 
mined  byDor  D'  and  H,  representing  pressures  "A"  and  "C",  and 
densities  J,  Lor  K.  As  points  D  or  D'  are  not  easily  obtained,  use  is 
made  of  points  F  or  G,  respectively.  In  the  case  of  curve  D'-H  the 
average  density  for  a  pressure  change  from  "C"  to  "A"  can  be  expressed 
as  follows; 


average  2  2  2  2 


The  average  density  for  curve  D-H  can  be  expressed  similarly. 


Substituting  Equation  (16)  into  Equation  (8) 


In  the  case  of  Fuel  A  (JP-4),  a  further  correction  is  required  due  to 
the  vapor  pressure  of  the  fuel.  Fuel  B  (KJ-  ll,  having  a  vapor 
pressure  below  0.  1  psia  at  ambient  temperature,  was  not  corrected 
for  this  factor. 

The  presence  of  appreciable  percentages  of  fuel  vapors  in  the  ullage 
requires  some  further  modification  of  the  basic  relationships  derived 
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in  Equation  (8)  to  take  into  account  the  partial  preeaure  of  the  fuel. 
Again  rigorous  calculations  could  not  be  used  due  to  practical  limffa' 
tions;  therefore  certain  approximations  were  used.  Dalton's  Law  of 
additive  pressure  is  one  of  the  approximations  assumed. 


Pf  +  Px 


The  fuel- vapor  correction  would  be  relatively  simple  to  make  if  the 
vapor  pressure  was  not  a  function  of  temperature  and  equilibrium 
achievement  time.  The  equilibrium  deviation  is  due  to  the  fact  that 
after  the  tank  is  vented  (as  in  a  climb  to  altitude  at  constant  guage 
pressure))  time  is  required  for  the  liquid  fuel  to  exert  its  full  vapor 
pressure.  This  is  similar  to  ordinary  evolution.  The  process  of  trans 
ferring  molecules  from  the  liquid  state  to  the  vapor  state  is  not  instan>a 
neous.  The  process  requires  a  finite  time. 

The  total  molal  quantities  entering  the  reservoir  are  approximately 
defined  by 


An^  = 

^"v.x*^"v.f>^"e 

(20) 

Ane  = 

(21) 

An  = 

^Pv.,V,  ^ 

e 

RTj.  Zy  RT^ 

(22) 

When  the  vapor  pressure  is  zero,  ^  is  zero.  If  the  vapor  pressure 
is  finite,  Anf^  y  increases  to  a  limit  as  a  function  of  time. 

Evaluation  of  Anf  y  is  involved.  Certain  simplifying  assumptions  can  be 
made.  The  two  extreme  cases,  assuming  both  to  be  without  gas  evolution 
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would  be  the  following  examples; 

Case  1  -  Vapor  pressure  constant  during  a  pressure  change  (climb). 

This  condition  would  be  closely  approximated  in  a  very 
slow  pressure  drop. 

Case  2  -  A  relatively  long  vapor-pressure-equilibrium  time 
causing  the  fuel-partial-pressure  to  drop  in  direct 
proportion  to  the  pressure  decrease.  This  condition 
would  be  approximated  in  very  rapid  pressure  drops. 

In  case  1,  the  ratio  of  moles  of  fuel  vapor  to  total  moles  m  the  reservoir 
would  equal  the  ratio  of  the  fuel-vapor-partial  -  pressure  to  the  total 
pressure.  (Ideal  conditions  are  appllrable»  inasmuch  as  the  reservoir 
pressure  is  subatmo spheric).  Further  the  molal  ratio  would  equal  the 
ratio  of  the  vapor  pressure  to  the  arithmetic  average  of  the  tank  pressure 
before  and  after  a  pressure  decrease. 


Pf.t 


(2  3) 


In  case  2,  the  ratio  of  moles  of  fuel  vapor  to  total  moles  m  the  reservoir 
is  equal  to  the  ratio  of  the  respective  pressures.  The  conditions  stated 
for  case  2  are  quite  similar  to  those  of  a  tank  containing  only  pressuriz- 
ii^  gas  and  fuel  vapor  with  little,  if  any,  liquid  fuel  present.  Therefore, 
the  quantity  of  gas  leaving  the  tank  is  directly  proportional  to  the  pressure 
change.  The  constituents  leave  the  tank  in  the  same  proportions  as  thev 
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exist  in  the  tank. 


AP. 


^t.  1  "  ^t.  2 


Ap £  ^  represents  the  partial  pressure  change  of  fuel  vapor  due  to  the 

change  in  total  pressure.  The  smaller  the  value  of  the  tank  pressure, 
the  larger  the  value  of  Ap  £  ^  ^  fixed  pressure  decrease.  The 

maximum  value  of  Ap^  ^  is  equal  to  the  vapor  pressure. 

For  Fuel  A  (JP>4)f  the  vapor  pressure  is  approximately  1.  4  psia  at 
ambient  conditions.  The  maximum  pressure  increment  used  was  IS 
psi  and  the  minimum  pressure  was  20  psia.  Under  these  conditions, 
the  maximum  value  of  Ap^  ^  can  be  determined  as  shown  below, 
assuming  case  2  conditions  exist. 


^f.  t,  1 
*x,  t,  1 


^f.t.2 

^x,  t,  2 


where 


j  =  1.4  P  xt  I  “  35  -  1.  4  =  33.  6  psia 


p  ,  _  =  Unknown  p  ^  =  20  - 

f  '  *^xt,  2 


Pft.2 


Solving  Equation  (26)  by  substitution  gives  a  value  of  p  ^  2  ~ 


=  Pf.tl  -'■f.U  = 


-bADD  TR  60-767 


62 


This  represents  the  maximum  deviation  of  vapor  pressure  than  can 
occur  in  the  test  system.  Since  the  pressure  change  is  linear,  an 
average  vapor  pressure  can  be  used  satisfactorily. 


An  -  Pt  Avg  p 

-f,  r  _  f,  r  _ 


In  the  example  given,  the  average  vapor  pressure  (p^  is 


1.4  -  0,  8 


=  1. 1  psia. 


Comparing  Equation  (23)  with  Equation  (28)  one  can  see  that  the 
deviation  in  vapor  pressure  and  thus  moles  of  fuel  vapor  entering  the 
reservoir  might  vary  by  a  factor  of  roughly  25%.  In  actuality,  the 
variation  is  much  less  since  neither  extreme  is  realized  in  practice. 
Further,  since  the  An.  is  a  small  fraction  of  n  ,  it  was  felt  that  the 

X I  f  IT 

variation  was  within  a  margin  of  experimental  error.  It  was  therefore 
arbitrarily  assumed  that  case  1  conditions  were  fully  applic  able  and 
the  relations  derived  therefrom  were  used.  Combining  Equation  i22» 
and  Equation  (24),  we  obtain 


AP  V  P  V 
A  r  r  V,  X  V 

^”e  ^  RT  '  Z  RT 
r  V  V 


2p,  AP  V 

^f,  t  r  r 

.  +  P  -  RT 
t,  1  t,  2  r 


X  =  Px,  1  -  Px.  2  =  <*■..  1  -  Pf'  -  'Pt.  2  -  Pf>  '^”1 


V  (30) 


An  = 
e 


AP  V 
r  r 


AP  V 

V  V 

Z  RT 

V  V 


^Pf,.  ~l 


TH  60-767 


63 


Combining  Equation  (16),  bringing  in  the  time  and  weight  factorfl 
results  in  the  overall  relationship 


The  compressibility  factor  (Z^)  is  the  mean  compressibility  of  the 
vapor  mixture.  was  assumed  to  be  equal  to  the  compressibility  of  the 
pressurizing  gas  as,  in  the  regions  where  it  becomes  appreciable,  the 
diluting  effect  of  the  fuel  vapor  is  negligible. 

CONCLUSION: 

For  ethane,  the  relationship  shown  in  Equation  (32)  was  used.  Equation 
(12)  was  used  for  the  air  and  nitrogen  in  Fuel  B.  The  air,  nitrogen- 
Fuel  A  system  required  a  further  modification  of  Equation  (32).  The 
expansion  coefficient  OC  was  used  in  place  of  densities. 
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Experimental  studies  were  performed.  The  apparatus  used  to  make  the 
studies  is  shown  in  Figures  26  through  29*  All  runs  were  made 
at  ambient  temperatures  (50"F  to  90®F).  Agitation  was  accomplished 
by  the  rocking  mechanism  whose  speed  could  be  varied  by  changing  a 
pulley  ratio.  The  agitation  amplitude  was  held  constant  at  approximate¬ 
ly  j  inch  as  measured  at  the  edge  of  the  tank  oven.  A  detailed  descrip¬ 
tion  of  the  equipment  and  the  procedure  used  is  given  in  the  paragraph.s 
entitled ’’Test  Equipment  and  Procedure'*.  Details  of  the  method  used 
for  determining  evolution  rate  arc  given  in  the  calculations.  Fijiures 
53  through  44  illustrate  the  experimental  re.sults. 

The  following  general  statements  can  be  made  about  the  air-fuel 
system: 

1.  Evolution  rate  is  a  linear  function  of  pressure-chaugc-r.ite 

2.  Increased  agitation  increases  ti  e  r.ite  of  evolution. 

The  effect  of  air-pressure-level  on  evolution  rate  does  not 
appear  to  be  appreciable,  at  least  to  t'O  psi. 

4.  Evolution  rate  is  a  function  of  . ijuilibrlum  solubility,  i.  e.  , 
the  greater  the  solubilifv  the  urcater  the  evolution  rate. 

5.  Deviation  from  equilibrium  increases  as  'He  ra'e -of-pr'  .  ui  ■  - 
change  increases. 

The  nitrogen-fuel  system  was  not  investigated  very  thoroughly.  It  is 
expected  that  the  generalized  statements  made  concerning  the  air  fuel 
system  are  applicable  to  the  nitrogen-fuel  system. 

The  ethane-fuel  system  differs  from  the  air-fuel  system  in  one  respe<  t 
the  ethane  pressure-level  greatly  effects  the  evolution  rate.  This  mav 
be  due,  in  part,  to  the  deviation  from  ideality  of  ethane.  Another 
factor  in  the  ethane-fuel  system  that  may  affect  the  evolution  rate  is 
the  rather  extreme  change  of  viscosity  and  presumably,  surface  tension 
with  solubility;  i.  e.  ,  with  partial  pressure. 


A  comparison  of  the  air-JP-4  and  the  ethane- JP-4  fuel-systems,  at  low 
rates  of  pressure  change^  shows  that  the  evolution  rate  of  air  is 
approximately  |  that  of  ethane.  The  solubility  of  air  is,  however,  1/20 
that  of  ethane  under  the  same  conditions.  This  infers,  of  course,  that 
air  can  be  removed  from  the  fuel  much  more  readily  than  ethane. 

The  effect  of  agitation  on  the  fuel  systems  was  fairly  well  defined 
except  for  the  case  of  zero  agitation.  Figure  39  shows  the  extreme 
variation  of  test  results.  It  is  believed  that  the  large  deviations  are 
due  to  our  inability  to  obtain  equilibrium  under  conditions  of  zero 
agitation. 

Although  all  runs  were  made  at  ambient  temperature,  the  temperature 
varied  approximately  5®F  to  lO'T*  during  a  test.  Weather  variations 
dui-ing  the  test  period  covered  the  temperature  range  of  SO^F  to  90*^, 
averaging  roughly  70°F.  These  temperature  variations  are  one  cause 
of  minor  deviations  noted  on  the  data  plots.  The  effect  of  temperature 
on  evolution  rate  is  quite  appreciable  in  the  ethane-fuel  system  due  to 
the  large  effect  of  temperature  on  ethane  solubility  and  thus,  fuel 
viscosity  and  surface  tension. 

The  air-fuel  system  results  obtained  in  this  test  are  not  in  agreement 
with  those  of  Schweitzer  and  Szebehely.  This  is  believed  due,  in  part, 
to  the  difference  in  agitation  severity  and  the  vapor-liquid-ratio  used. 
Schweitzer  and  Szebehely  used  a  shake-table  type  agitator  with  a  1-inch 
stroke  and  a  vapor-liquid- ratio  of  4:1.  Work  done  on  the  program 
reported  herein  was  accomplished  with  a  rocking  agitator  with  a  j-inch 
maximum  stroke  and  a  vapor  -  liquid-  ratio  of  approximately  l  iO 
Further  information  is  given  in  Appendix  I  of  this  Part,  on  the 
work  of  Schweitzer  et.  al.  The  equation^  of  Appendix  I  were  derived 
by  Staff  Members  of  the  Stanford  Research  Institute  as  part  of  a  sub¬ 
contract  under  the  literature  survey  portion  of  this  program 

The  ethane-fuel  system  is  also  not  in  agreement  with  the  results 
predicted  by  Schweitzer  et  al  In  this  rase  an  additional  factor  adds 
to  the  discrepancy;  ethane  does  not  follow  Henry's  Law  as  is  required 
in  the  calculations  of  Schweitzer  and  Zebehely  This  factor  shows  up 
in  the  effect  of  pressure  level  on  evolution  rate 
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APPENDIX  I  DISSOLUTION  AND  EVOLUTION  OF  GASES  IN  LIQUIDS 
(Developed  by  Physical  Chemistry  Staff  of 
Stanford  Research  Institute) 


The  following  equations  are  based  on  the  work  of  Schweitzer,  et  al. ,  >}' 
who  found  that  the  rate  of  evolution  or  dissolution  of  a  gas  from  various 
hydrocarbon  fractions  was  proportional  to  the  degree  of  supersaturation 


Assume 


dv 
dt  "" 


-k  (v-v^) 


(1) 


where  K=  pseudo  first  order  rate  constant,  which  is  a 
function  of  the  direction  of  the  process 
(except  very  close  to  equilibrium),  temperature, 
nature  of  the  gas  and  of  the  liquid,  rate  of 
agitation,  etc.  ,  but  not  of  the  pressure. 


♦  Schweitzer,  P.  H.  and  Szebehely,  V.  G.  Gas  Evolution  in 

Liquids  and  Cavitation.  Journal  of  .Applied  Physics.  Volume  21. 
December,  1950.  pp.  1218  -  12  24. 
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V 


volume  of  gas  dissolved  (measured  at  NTP)  per 
volume  of  liquid,  at  time  t, 

Vg  =  volume  of  gas  dissolved  per  volume  of  liquid 
that  would  be  in  equilibrium  with  the  pressure 
p  at  time  t. 

Also  assume  that  the  dissolved  gas  is  sufficiently  dilute  so  that 
Henry's  law  is  obeyed: 

Ve  =  hp  (2) 

Case  I.  Closed  system,  constant  temperature.  At  t  =  o,  v  =  Vq 
and  p  =  Pg  (not  in  equil.  with  Vq). 

The  material  balance  is  given  by 


pVg  +  vVi  =  pf  Vg  +  V£  Vi 


p  -  Pf  =  r(vf  -  v) 


(3) 


where  p£ 
vf  = 


Vl 


final  equil.  pressure  at  t  =  '-O 

final  equil,  volume  /  volume  in  the  liquid 

total  volume  of  gas  phase 

total  volume  of  liquid  phase 

volume  of  liquid  /  volume  of  gas 


Substituting  (2)  and  (3)  into  (1)  (and  noting  V£  =  hp£) 


dv 

3  -kv  +  kvQ  +  k  cCt 


(8) 


Equation  (8)  ia  a  linear  equation  of  1st  order  with  constant  coefficients 
whose  solution  is 


s  Ce"*'*  +c<Lt  -  ^  (koC-  Vq) 


=  Ge  ***  +  ©Ct  +  X, 


At  t  =  o,  V  =  Vg, 


...  C  =  -2^  and 

k 


v-v^soCt  +  SL  -1) 


(9) 


Noting  that  Vq  -  =  -oCt  (eq.  (7)),  equation  (9)  can  be  put  in  the 

somewhat  more  useful  form: 


_2 -  =  I  -  j^(  1  -  e-*'*) 

v  -V.  kt 


(9a) 


The  quantity  (vq  -  v)  /  (Vq  -  Vg)  is  the  ratio  of  the  gas  actually  evolvec 
(or  absorbed)  to  that  which  would  be  evolved  at  time  t  had  equilibrium 
obtained.  Equation  (9)  can  also  be  written  in  terms  of  pressure: 


Pe  -P 


(9b) 
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H  one  interested  in,  or  forced  to  measure,  the  total  loss  of  gas  from 
an  open  system  (e.  g. ,  the  loss  from  a  tank  partially  filled  with 
gasohne  and  equipped  with  a  vent)  the  equation  must  be  recast  slfghtiy 
Denoting  the  gas  lost  (measured  at  NTP)  by  G,  a  material  halanr  e  give-. 


-iSf-  +  V, 
dt  ^  dt 


=  -VgiE- 

»  dt 


=  -v,  ^  -  V.  Je 

dt  ^  dt  g  dt 


Sub.tituting  erom  equation  (8)  and  noting  -hz  and  retailing a  -a. 


=  kVj  (V  -  Vo  -^t)  -  ^ 


=  °<.Vi(e-^*  -  1) - S. 


Integrating 


V,e' 


V  eC  t 


-  V,t  - 


c<V 

At  t  a  O,  G  =  O,  .  C  =  _ L 

k 


C<V  , 

80  G  =  _ 1, 

k 


(1  -  e-l  -  'jVZli 

h 


In  a  somewhat  more  convenient  form: 


=  1  - 


.  1 


Vg  f  V^h 


•(1  -  e'^t) 


(10a) 


where  Gg  is  the  total  loss  of  gas  if  equilibrium  had 
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SYMBOLS 


E 

M 

Evolution  rate 

P 

m 

Partial  pressure 

P 

m 

Pressure 

R 

M 

Gas  constant 

t 

3 

Time 

T 

= 

Temperature 

V 

3 

Volume 

W 

3 

Weight 

Z 

3 

Compressibility 

tX, 

3 

Liquid  expansion  factor 

A 

3 

Change- 

3 

Density 
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SUBSCRIPTS  AND  SUPERSCRIPTS 

Pressurizing  gas  evolving  from  liquid  fuel 

Fuel  vapor 

Initial  liquid 

Liquid  fuel 

Reservoir 

System  less  reservoir 
Tank 

Tank  vapor  space 
Pi’ossurizing  gas 
Change 
Condition  1 
Condition  2 


Condition  3 


26  «  Sehematle  of  Test  Tank  and  Accessories 


of  Gas 
Second 


Avcra.gt  iUt«  of  C*i  Evolution- -lb.  mole«  ofgai| 

lb.  fuel  ■econd 

Figure  36  •  Air  •  JP-4  Evolution  Rate  (150  CPM) 


Average  Rate  of  Gas  Evolution- -lb.  moles  of  gae 

_  lb,  fuel  second 

ilgure  37  -  Air  -  JP  -  4  Evolution  Rate  (240  CPM) 
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Average  Rate  of  Gas  Evolution--Ib.  moles  of  gas 

lb.  fuel  second 

Figure  39  -  Air  -  RJ-1  Evolution  Rate  (0  CPM) 


Average  Rate  oi  Gae  Evolution- -lb.  moles  of  ga» 

lb.  ioel  second 

Figure  i^O  -  Air  -  RJ-1  Evolution  Rate  (150  CPM) 


Average  Rate  of  Gas  Evolutiln-tlb.  Ul^a  of  f.as  ' 

lb.  fuel  second 

kl  -  Air  -  Rj-1  Evolution  Rate  (240  CPM) 


PART  III 


VAPORIZATION  CYCLE  AND  THERMODYNAMIC 
PROPERTIES  OF  JP-4,  RJ-1  AND  A 
SPECIAL  FUEL  DESIGNATED  JP-6-H 


INTRODUCTION 

This  report  was  first  issued  in  January,  1958  as  one  of  a  series  of 
Interim  Reports.  The  interim  Reports  were  submitted  to  Wright  Air 
Development  Center  by  Thompson  Products,  Inc.  West  Coast  Laboratory 
in  partial  fulfillment  of  research  program  No.  AF  33(616)3729. 

The  material  reported  here  covers  the  vaporization  cycle  and 
thermodynamic  properties  of  JP-4,  RJ-1  and  a  special  fuel  designated 
JP-6-H.  This  data  represents  results  of  completed  studies. 


OBJECTIVE 

Studies  were  made  under  this  part  of  the  program  for  two  basic  reasons: 

1)  The  determination  of  vaporization  characteristics  of  a 
JP-4,  a  special  JP-6,  and  an  RJ-1*  fuel  as  a  function  of 
pressure  level.  If  possible,  the  pressure  range  of  0  to  200 
psia  was  to  be  covered.  It  was  also  desired  to  determine 
the  specific  gravity  and  molecular  weight  of  these  fuels  as 
a  function  of  per  cent  vaporized,  and  the  variation  of  these 
physical  properties  with  pressure  level. 

2)  The  determination  of  enthalpy,  entropy,  specific  heat,  and 
heat  of  vaporization  of  the  three  fuels  listed,  as  a  function 
of  temperature  and  pressure.  As  much  of  this  information 
as  possible  was  to  be  presented  in  the  form  of  enthalpy- 
entropy  diagrams. 


*  Formerly  designated  Shell  UMF,  Grade  C 


SUMMARY 

Objectives  1  and  2  were  studied,  A  brief  summary  of  the  results 
achieved  and  the  determinations  made  is  given  in  the  following  para¬ 
graphs. 

Objective  1 

The  experimental  vaporization  work  was  accomplished  using  equipment 
(See  Figure  45  )  designed  to  obtain  a  minimum  of  refluxing  such  as 

would  be  experienced  in  aircraft  fuel  tanks  in  which  fuel  was  being  vapor¬ 
ized.  Results  are  shown  in  Figures  4b  through  S4*  figures  include 

calculated  equilibrium  vaporization  curves  and  ASTM  distillation  curves 
for  comparison  with  the  minimum  reflux  vaporization  data. 

Objective  2 

Experimental  specific  heat  results  are  shown  in  Figures  33  through  37* 

As  expected,  the  liquid  specific  heat  follows  the  C  =  A  +  BT  type  relation 
where:  C  represents  specific  heat,  A  and  B  are  constants,  and  T 
represents  temperature. 

The  enthalpy,  entropy  and  heats  of  vaporization  of  each  of  three  fuels 
studied,  were  computed  using  empirical  correlations  plus  experimental 
distillation  and  density  data.  Enthalpy  -  temperature  plots  in  the  O'*  to 
900°F  range  with  pressure  as  a  parameter  are  shown  in  Figures  58 
through  6o  Temperature  -  entropy  plots  in  the  same  temperature  range 
with  enthalpy  and  pressure  parameters  are  shown  in  Figures  fal  tSrcugh  63 

METHODS,  EQUIPMENT  AND  PROCEDURES 

Minimum  reflux  vaporizations  of  the  three  test  fuels  were  accomplished 
in  the  test  cell  shown  in  Figure  45*  The  fuel  vaporized  was  contained 
in  the  inner  tube.  A  heating  medium  (Ondina  oil)  was  circulated  through 
the  annular  space  between  the  tubes.  This  method  of  heating  was  chosen 
in  order  to  maintain  the  inner  tube  close  to  an  isothermal  condition. 
Temperature  stabilization  prevented  condensation  of  vaporized  constitu¬ 
ents  and  kept  refluxing  to  a  minimum.  The  only  area  where  refluxing 
was  possible  was  in  the  exit  tube.  A  capillary  tube  of  minimum  length 
was  used  to  minimize  the  effects  of  thi.s  refluxing. 
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Pressure  was  manually  controlled  by  a  valve  in  the  exit  line  and  was 
monitored  by  means  of  a  pressure  transducer.  The  temperature  was 
monitored  using  thermocouples.  Condensation  of  the  vaporized  products 
was  accomplished  by  means  of  a  wet  ice  bath.  Constant  evaporation 
rates  were  not  maintained  since  this  factor  does  not  affect  the  separation 
efficiency  of  the  test  apparatus,  providing  the  evaporation  rate  is  not 
sufficiently  rapid  to  cause  liquid  entrainment.  Several  tests  were  con¬ 
ducted  to  verify  separation  efficiency. 

Fractions  were  collected  in  ten  per-cent-by-volume  increments  at  each 
of  the  pressure  levels.  The  initial  and  final  temperatures  were  noted. 
The  density  of  the  collected  fractions  was  measured  using  a  calibrated 
pynenometer  and  a  precision  analytical  balance. 

Molecular  weight  determinations  were  made  for  each  of  the  fractions. 
Freeze^point  depression  methods,  with  benzene  as  a  solvent,  were 
used  to  make  these  determinations. 

A  Qomparison  Calorimeter,  shown  schematically  in  Figure  64  was 
used  to  determine  the  liquid  specific  heat  of  the  three  test  fueis. 

For  detailed  information  concerning  the  use  of  this  calorimeter  see 
Reference  No.  1.  In  addition  to  the  physical  analyses  and  ASTM 
distillations,  vapor  pressure  vs  temperature,  and  density  vs  temperature 
teats  were  made  on  the  fuels.  Specific  volume  techniques  and  the 
Specific  Volume  Apparatus  shown  in  Figire  65  were  used  to  make 
the  density  vs  temperature  determinations. 

The  Specific  Volume  Apparatus,  slightly  modified,  was  used  to  measure 
vapor  pressures.  A  vacuum  pump  and  the  liquid  nitrogen  condenser  in 
the  Specific  Volume  Apparatus  were  used  to  de-aerate  the  fuels.  After 
fuel  de-aeration,  the  vapor  pressure  was  measured  at  various  points  in 
the  50'F  to  SOO'F  temperature  range  by  means  of  ^  pressure  transducer. 
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ENTROPY  AND  ENTHALPY  CALCULATIONS 


The  Entropy  and  Enthalpy  charts  presented  in  Figure  ^  through 
Figure  63  were  calculated  in  the  following  manner: 

From  the  ASTM  plots  of  the  fuels  tested,  the  equilibrium  vaporization 
curves  were  obtained  (method  of  Edmister  and  Pollock,  see  Reference 
2).  The  information  obtained  was  extrapolated  to  low  pressures. 
Liquid  enthalpies  were  computed  using  basic  thermodynamic  relation¬ 
ships  and  Cragoe's  empirical  relation  for  specific  heat  (Reference  3). 


dH  = 

CpdT 

(1) 

II 

Hq  ^  y  c  dT 

Tj  P 

(2) 

C,  = 

0.  388  +  0.  00045T 

(3) 

'-L 

^0.  5 

Assuming  a  datum  temperature  of  0“F  i.  e.  ,  H"  =  O  at  0"F,  substituting^ 
(3)  and  integrating  gives 

2 

„  ^  0.  388T  +  0.  000255T 

^0.5  (4) 
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Evaluation  of  Equation  (4)  at  various  temperature  levels  results  in 
enthalpy  values.  These  values  when  plotted,  result  in  points  from  which 
the  saturated  liquid  line  of  an  enthalpy  temperature  chart  may  be  drawn. 
Ibobtain  the  saturated  vapor  line  and  complete  the  phase  envelope,  heats 
of  vaporization  must  be  determined. 

Heat  of  vaporization  is  defined  as  the  enthalpy  difference  between  the 
saturated  vapor  and  saturated  liquid  at  constant  temperature.  For  the 
fuels  studied  here,  this  constant  temperature  is  assumed  to  be  the 
mean  average  boiling;  point.  Correlations  such  as  those  shown  in 
Appendix  A  (Part  V,  Supplement  V  of  this  Report)  were  used  to  obtain 
the  heats  of  vaporization  at  one  atmosphere  and  at  the  mean  average 
boiling  points  Heats  of  vaporization  at  other  temperatures  were  obtained 
using  the  correlation  of  Watson(see  Reference  4). 


The 


A 


function  is  the  ratio  of  the  heat  af  vaporization  at  the 


desired  temperature  to  the  heat  of  vaporization  at  the  mean  average 
boiling  point.  The  addition  of  the  heat  of  vaporization  to  the  enthalpy 
of  the  saturated  liquid  at  the  mean  temperature,  results  in  the  enthalpy 
of  the  saturated  vapor. 


H 


sv 
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A  timilar  procedure  is  used  to  calculate  the  per  cent  vaporized  lines. 
Since  the  enthalpy  is  normally  plotted  on  a  weight  basis  and  the  per  cent 
vaporized  on  the  equilibrium  vaporization  curve  on  a  volume  basis, 
suitable  corrections  must  be  made.  To  make  the  corrections  for  Fuel  A, 
density  measurements  from  experimental  minimum-reflux  vaporizations 
were  used.  Knowing  the  mean  temperature  of  the  fuel  at  different 
percentages  of  vaporization  and  the  enthalpy  envelope,  allows  the  evalua¬ 
tion  of  enthalpy  change  for  each  of  the  percentages  vaporized.  In  this 
statement  the  mean  temperature  is  presumed  to  be  the  arithmetic  average 
boiling  point  of  the  fraction,  i.  e. , 


*^mean  0-20  = 


^20  »  "O 
2 


'■mean  20-40  s 


"^40  *  ^20 


",  etc. 


In  the  case  of  narrow  boiling  range  fuels,  such  as  fuels  B  and  C,  these 
corrections  are  negligible. 

At  the  mean  temperature,  the  heat  of  vaporization  is  obtained  for  the  one- 
atmo^phere  curve  using  the  enthalpy  differential  between  the  saturated 
vapor  and  saturated  liquid  line.  If  the  enthalpy  of  the  twenty  per  cent 
fraction  is  being  computed,  twenty  per  cent  of  the  enthalpy  differential 
is  added  to  the  liquid  enthalpy  at  the  end  point  temperature;  at  forty  per 
cent,  forty  per  cent  of  the  differential  is  added,  etc.  At  another  pressure 
level  the  assumption  is  made  that  the  ratio  of  the  enthalpy  increments, 
for  any  fraction,  will  change  as  the  ratio  of  the  total  enthalpies  change. 


@  14. 7  _  ^^ZO-40  g  14.7 

(§  p  AH20-40  @  p 


If  the  region  being  investigated  is  not  near  critical  conditions  this  is  a 
reasonable  assumption.  The  enthalpy  of  the  saturated  vapor  is  equal  to: 
1}  the  summation  of  the  liquid  enthalpy  (at  the  initial  boiling  point),  and 
2)  the  enthalpy  of  vaporization  of  all  increments. 
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(7) 


To  evaluate  the  specific  heat  at  constant  pressure  of  the  vapors  in  the 
superheat  region,  the  residual  method  of  Deming  and  Shupe(References 
5  and  6),  plus  generalized  compressibility  and  specific  heat  charts, 
were  used. 


ZRT 

P 


(13) 
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The  second  differential 


ial 


3^  V 
3  T' 


was  evaluated  from  Equation  (14) 


by  double  graphical  differentiation  at  various  pressure  levels.  Graphical 
integration  of  the  results  and  substitution  into  Equation  (11)  gives  the 
correction  to  the  specific  heat  due  to  pressure. 

As  before,  the  enthalpy  change  due  to  increasing  temperature  is  the 
integral  of  the  specific  heat  equation.  In  the  case  of  the  vapors,  to 
simplify  calculations,  use  of  the  mean  specific  heat  was  made. 


=  C  (AT) 

p  mean 


The  mean  specific  heat  was  taken  as  the  arithmetic  average  over  a 
small  temperature  increment.  The  AH  value  computed  then  was 
added  to  the  saturated  vapor  enthalpy  and  the  enthalpy-temperature 
diagram  was  completed. 

Entropy  computations  were  based  also  on  a  datum  of  0“F.  The 
thermodynamic  relationships  used  were  : 


C  dT 
P 
T 
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To  compute  the  saturated  liquid  line  on  the  temperature-entropy  chart, 
the  Cragoe  equation  was  substituted  in  Equation  (16),  The  equation 
was  then  integrated,  resulting  in 


181  (In  T2 


-  In  460)+  0.  00045C 


(18) 


This  equation  represents  the  saturated  liquid  line. 

The  temperature-entropy  charts  presented,  (see  Figures  through 

63  were  derived  from  Equations(l6)  and  (17)  using  the  values  of 
change-in-enthalpy  presented  in  the  temperature-enthalpy  charts.  See 
Figures  38  through  60« 

In  the  superheat  region  each  temperature,  substituted  in  Equation  (16), 
was  the  arithmetic  average  which  closely  approximates  the  mean  over 
the  small  increments  taken. 


Critical  constants  were  determined  using  correlations  contained  in 
Part  V,  Supplement  V.  Data  in  the  vicinity  of  the  critical  region  were 
obtained  by  ext  rapolation. 

DISCUSSION 


A  study  of  the  vaporization  cycle  and  thermodynamic  properties  of 
three  test  fuels,  a  JP-4,  a  special  JP-6-H  and  an  RJ-1=^'  was  completed. 
The  research  study  results  are  contained  in  the  data  presented  in  this 
report  and  concludes  the  investigation. 

The  RJ-1  fuel  tested  is  in  effect  a  Shell  UMF  Grade  C  and  conforms  to 
specification  MIL-F-25558.  Both  the  JP-4  and  th^  RJ-1  test  fuels 
appeared  to  be  representative  samples.  The  JP-6-H  test  fuel  mee’s 
specification  MIL-F-25656.  This  fuel  is  considered  special  because 
of  the  very  narrow  boiling  range.  Typical  JP-6  fuel  boils  from  280°F 
to  530*^.  This  JP-6-H  boiled  from  328®F  to  405°F.  The  physical 
properties  of  the  test  fuels  are  discussed  in  the  following  paragraphs. 

♦  Formerly  designated  Shell  UMF  Grade  C. 
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Test  Fuel  A  (JP-4)  and  Test  Fuel  B(RJ-l)  were  purchased  from  the 
Shell  Oil  Company.  Test  Fuel  C  (JP-6-H)  was  provided  by  Wright 
Air  Development  Center,  The  ASTM  distillation  plot  of  the  JP-4  fuel 
shows  a  rather  steep  rise  in  the  initial  ten>per-cent-distilled  range. 

Figures  66  through  68  show  the  variation  of  specific  gravity  with 
temperature.  For  fuels  A,  and  C  the  decrease  of  density  with 
increasing  temperature  is  nearly  linear  until  approximately  300°F. 

Fuel  B  shows  a  linear  change  of  density  over  the  entire  test  tempera¬ 
ture  range  (0^  to  500'*F). 

The  liquid  specific  heat  plots,  shown  in  Figures  55  through  57 
follow  the  normal  trend  for  hydrocarbon  fuels.  Fuel  A  follows  the 
equation  C  =  0.  472  f  0.  0005T;  Fuel  B  follows  the  equation  C  =  0,  420 
+  0.  00044T;  and  Fuel  C  follows  the  equation  C  =  0.  468  f  0.  00049T. 

Each  of  the  fuels  demonstrates  the  normal  petroleum- vapor-pressure 
characteristics,  i.  e. ,  plotting  vapor  pressure  versus  temperature  on.  a 
Cox  chart  results  in  a  straight  line.  See  Figures  69,  70  and  71. 

Minimum  reflux  distillation  data  at  14.  7  psia  (see  Figures  46-48 
when  compared  with  ASTM-distillation  data  (Figures  72-74  shows 
the  effect  of  reflux  on  the  vapor  product.  Minimum  reflux  distillations 
have  a  considerably  higher  initial  boiling  point  and  a  somewhat  loweV 
end  point  than  ASTM  distillations.  This  is  due  to  the  effect  of  reflux. 

In  aircraft  flying  at  speeds  sufficient  to  cause  boiling  in  the  fuel  tanks 
due  to  aerodynamic  heating,  the  conditions  in  the  tanks  will  closely 
follow  those  of  a  minimum  reflux  distillation.  This  statement  is 
predicated  on  the  assumption  that  the  tank  is  held  at  a  constant  pressure. 

Figures  49  through  54  sh»v  the  variation  of  density  and  mole¬ 
cular  weight  as  a  function  of  the  fuel  fraction  vaporized  for  minimum 
reflux  distillations.  The  fuel  fraction  vaporized  is  plotted  as  volume 
per  cent.  The  variation  between  volume  and  weight  per  cent  is  neg¬ 
ligible  when  the  distillation  range  is  200*^  or  less.  The  plots  of 
molecular  weight  and  density  flatten  out  with  increasing  pressure,  the 
e.xtremo  occurring  at  the  critical  pressure. 

The  equilibrium  vaporization  curves  have  been  included  for  additional 
information.  See  Figures  75  through  77»  These  curves  show 
the  vaporization  conditions  in  a  tank  where  pressure  has  suddenly  been 
lowered  without  the  loss  of  fv.el  vapor  froni  'he  system. 
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The  vapor  epccific  heat  data  shown  fn  Figures  78  through  80 
were  computed  from  generalized  Pre3s>:re- Volume-Temperature 
(P-V-T)  data  and  vapor  specific  heat  correlations  at  one-atmosphere, 
using  the  residual  method  of  Deming  and  Shupe.  See  Reference  5.  It 
was  assumed  that  deviations  due  to  pressure  helow  14.  7  psia  were 
insignificant.  Therefore,  the  specific  heat  at  14.  7  psia  was  used  for 
Cp“in  the  calculations. 

The  temperature-enthalpy  and  temperature-entropy  diagrams  (Figures 
58  through  63  )  were  computed  from  ASTM  distillation  data  and 

numerous  correlations.  The  calculations  involved  are  described  in 
this  Supplement  in  paragraphs  entitled  'Entropy  and  Enthalpy  Calcula¬ 
tions'.  These  plots  are  most  applicable  to  a  closed  fuel  system  from 
which  liquid  fuel  only  is  supplied  to  the  engine  and  where  the  fuel  tank 
vapors  are  not  vented.  In  the  case  of  a  wide  cut  distillate  Where  the 
tank  is  vented  or  vapors  are  supplied  to  the  engine,  the  usefulness  of 
the  charts  is  limited.  The  degree  of  limitation  depends  upon  the  exlent 
of  vapor  removal  from  the  fuel  tank.  A  number  of  enthalpy-temperature 
and  entropy-temperature  curves  can  exist  for  the  same  initial  fuel.  The 
manner  in  which  the  fuel  is  vaporized  will  determine  the  curve.  It  is 
pointed  out  that  changes  in  fuel  composition  due  to  distillation  vary  with 
pressure  level.  (See  Figures  49  through  34  )-  Therefore  for 

ultimate  use,  three  dimensional  charts  (with  the  change  in  composition 
of  the  hydrocarbons  added)  should  be  constructed  for  specific  applica¬ 
tions.  In  lieu  of  a  three  dimensional  chart,  a  series  of  charts  covering 
successive  fuel  residues  is  required  to  define  enthalpy  and  entropy 
changes. 

CONCLUSION 

Although  the  information  shown  on  the  various  plots  included  in  this 
report  is  valid  only  for  the  three  fuels  tested,  certain  generalizations 
can  be  made.  In  the  case  of  Fuel  A,  a  representative  JP-4  fuel,  use  of 
the  data  for  estimation  purposes  in  the  design  of  systems  and  equipment 
should  be  satisfactory.  In  the  case  of  Fusl  B,  the  deviation  between 
fuels  of  this  type  being  minor,  use  of  the  data  for  most  cases  should  be 
well  within  normal  engineering  tolerances.  However,  the  data  for  Fuel 
C  should  be  treated  with  caution  as  it  is  not  a  typical  JP-6.  Fuel  C  is 
classified  as  a  'heart  cut',  i.  e. ,  a  narroM’  boiling  range  fuel. 
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SYMBOLS 


C  »  Specific  Heat  BTU/lb 

H  *  Enthalpy  -  BTU/lb 

J  s  Mechanical  equivalent  of  heat 

K  s  Characterization  factor 

MW  s  Molecular  Weight 

P  »  Pressure  -  psia 

Q  s  Heat  of  Combustion  BTU/lb 

R  s  Gas  constant  -  10.73  psia,  ft^/*^,  lb-mole 

S  =  Entropy  -  BTU/lb 

T  =  Temperature  -  ^  or 

V  a  Volume  -  ft^ 

Z  a  Compressibility  factor 
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CC  =  Volume  residual 
A  =  Increment 

X  =  Heat  of  vaporization  -  BTU/lb* 

=  Density  at  60‘*F 

SUBSCRIPTS  AND  SUPERSCRIPTS 
c  =  Critical 

k  s  Normal  mean  average  temperature 

L  -  Saturated  liquid 

m  =  Mean 

p  =  Constant  pressure 

r  =  Reduced  conditions 

sv  =  Saturated  vapor 

t  =  Total 

V  =  Vaporization 

0  to  20,  20  to  40,  etc.  =  0-20%,  20-40%  Fractions 

1  =  Initial  condition 

2  =  Final  condition 

*  Also  used  to  indicate  heat  of  vaporization. 
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FIGURE  66.  SPECIFIC  GRAVIT Y  -  FUEL  A  ,  J P -  4 
WADD  TR  60-767  129 


SPECIFIC  GRAVITY 


SPECIFIC  GRAVITY  -  TEMPERATURE  PLOT 


FIGURE  68.  SPECIFIC  GRAVITY,  FUELC,  JP-6-  H 
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FIGURE  72.  DISTILLATION  PLOT,  FUEL  A,  JP 
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FIGURE  73.  DISTILLATION  PLOT,  FUEL  B 
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PART  IV 


VAPORIZATION  CYCLE  AND  THERMODYNAMIC 
PROPERTIES  OF  DECALIN 


INTRODUCTION 

A  report,  Ing.  Er.  232  ,  a  study  of  the  thermodynamic  properties 
of  decalin  was  fourth  in  a  series  of  Interim  Reports.  Each  of  the  Interim 
Reports  was  offered  in  partial  fulfillment  of  a  research  contract 
conducted  under  United  States  Air  Force  Contract  AF  33(616)3729. 

The  results  contained  in  Part  IV  of  this  report  Ing.  Er.  234,  are 
concerned  with  the  completed  study  of  the  thermodynamic  properties  of 
decalin. 

OBJECTIVE 

The  purpose  of  this  phase  of  the  program  was  to  determine  the  thermo¬ 
dynamic  properties  of  commercial  dec  ilin.  Emphasis  was  placed  on 
the  preparation  of  enthalpy-temperature  and  entropy-temperature 
diagrams.  In  addition  to  the  thermodynamic  properties  such  as  specific 
heat,  vapor  pressure,  and  ASTM  distillation;  the  effect  of  pressure  on 
the  applicable  properties  was  to  be  determined. 


SUMMARY 


Decalin  ,  or  more  precisely,  decahydronapthalene,  is  a  hydrocarbon. 
Decalin  possesses  mahy  characteristics  which  make  it  excellent  for 
use  as  an  aircraft  fuel,  at  elevated  temperatures.  It  has  a  relatively 
high  boiling  point  and  should  be  thermally  stable  to  temperatures 
approaching  800°F.  The  commercially  available  product  consists,  in 
the  main,  of  a  mixture  of  the  cis  and  trans  isomers  along  with  a  small 
quantity  of  impurities.  The  boiling  range  is  appro-xima tely  30°F  as 
determined  by  an  ASTM  distillation.  (This  represents  a  liquid  temper¬ 
ature  range  of  approximately  10‘F.) 
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Some  information  is  available  in  the  literature  References  1  through 
7.  Where  feasible,  this  information  was  used  to  corroborate  and 
complement  experimental  work  done  by  the  West  Qoast  Laboratory  of 
Thompson  Products,  Inc. 


The  thermodynamic  charts  presented  In  this  report  were  derived  in 
the  manner  described  in  Supplement  III  of  this  report,  with  a  few 
minor  modifications. 


DISCUSSION 


The  teat  fuel,  commercial  decalin,  was  supplied  by  Wright  Air 
Development  Center  (WADC).  Figure  81  shows  the  distillation  plot 
of  the  fuel  as  determined  by  ASTM  method  D- 158-54.  It  should  be 
realized  that  in  the  ASTM  method,  vapor  temperatures  are  measured. 
This  being  the  case,  superheating  of  the  vapors  occurs  quite  easily, 
giving  rise  to  the  comparatively  wide  variation  in  boiling  range.  Although 
a  distillation  range  of  32°F  is  shown  on  Figure  81  the  boiling  range  using 
liquid  temperatures  as  the  criterion,  shows  a  range  of  approximately 
10“F. 

The  effect  of  temperature!  on  specific  gravity  is  shown  in  Figure  82* 

Up  to  approximately  400®F,  the  density  decrease  is  linear  with  increasing 
temperature.  '  ■ 

A  plot  of  specific  heat  is  shown  in  Figure  63*  This  plot  is  typical  of 
the  hydrocarbons  and  can  be  described  mathematically  by  the  expression 

C  =  0.  351  +  6.  2  X  10-*  T  where  T  is  'T*  and  C  is  BTU/lb®F. 

P  P 

Flash  vaporization  curves.  See  Figure  8lv  were  determined  using 
the  method  of  Edmister  and  Pollock  (Reference  9).  It  is  readily  evident 
that  for  all  practical  purposes,  the  test  fuel  can  be  considered  pure. 

Figure  85  is  a  vapor  pressure  plot.  It  should  be  noted  that  it  also 
represents  the  bubble  point  curve  on  an  Equilibrium  Flash  Vaporization 
(E.  F.  V.  )  plot. 
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The  variation  of  vapor  specific  heat  with  temperature  and  pressure  is 
illustrated  in  Figure  66*  The  calculations  were  made  using  the 
method  described  in  Supplement  111  of  this  report  with  the  exception  that 
C  *or  the  specific  heat  at  one-atmosphere  was  evaluated  from  the 
enthalpy-temperature  diagram  shown  in  Figure  87*  (See  the 
Calculations  of  this  Supplement  for  additional  information). 

The  thermodynamic  charts.  Figures  87  and  88  were  determined 
by  methods  described  in  detail  in  fart  III  of  this  report#  These 
methods  are  abridged  in  the  calculation  section  of  Part  IV. 

Minimum  reflux  distillations  were  not  performed  on  decalin  due  to  the 
narrow  boiling  range.  The  E.  F.  V.  plot  shown  in  Figure  85  reflects 
quite  closely  what  would  occur  during  a  minimum  reflux  distillation  of 
a  fuel  with  a  narrow  boiling  range,  such  as  decalin. 

For  the  same  reason  i,e. ,  narrow  boiling  range,  the  variation  of 
condensed  vapor  density  and  molecular  weight  as  functions  of  per  cent 
distilled  and  pressure  were  not  experimentally  evaluated.  The  maximum 
variations  in  molecular  weight  would  in  all  probability  be  1 1%  which  is 
roughly  the  limit  of  our  present  test  equipment..  Condensed  vapor 
density  would  also  show  a  maximum  variation  of  ±  1%. 

Figure  89  shows  the  heat  of  vaporization.  The  apparatus  shown 
schematically  in  Figure  90  was  used  to  determine  the  value  at  one 
atmosphere.  The  relationship  of  K.  M.  Watson  (Reference  8)  was  used 
to  evaluate  the  heat  of  vaporization  over  the  range  shown. 
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CALCULATIONS 


The  mathematics  involved  in  the  determination  of  the  thermodynamic 
properties  of  decalin  are  summarized  as  follows: 


dH  =  C  dT 

n 
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(5) 


l£  further  information  is  desired,  the  reader  is  referred  to  Part 
in  of  this  report. 

For  decalin  the  C  ,  value  was  determined  from  the  relationship 
p«  L 

s  0.  3Bl  +  6..2  X  10"“*  T 

P 

The  value  of  C  *  was  determined  by  using  equation  (2)  for  the  enthalpy 
P 

of  the  saturated  vapor  line  in  the  pressure  range  below  14.  7  psia  and 
a  C  correlation  given  in  Part  V  of  this  report* 


All  the  terms  on  the  right  hand  side  of  equation  (6)  can  be  determined 

from  the  H-T  chart.  The  C  value  thus  determined,  in  the  low 

P-  g 

temperature,  was  checked  against  the  uncorrected  values  given  in  the 
correlation.  A  suitable  correction  factor  then  was  determined.  This 
method  was  utilized  in  place  of  the  UOP  K  correction  factor  which  is 
not  as  suitable  in  an  application  concerned  with  pure  compounds. 

Heats  of  vaporization  were  determined  both  experimentally  and  by  the 
use  of  correlations  where  data  was  not  available. 
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S  YMBOLS 


C  s  specific  heat 

H  =  enthalpy 

J  s  mechanical  equivalent  of  heat 

K  =  characterization  factor 

P  =  pressure 

S  =  entropy 

T  =  thermodynamic  temperature 

V  =  volume 

SUBSCRIPTS  AND  SUPERSCRIPTS 

L  -  saturated  liquid 

g  =  vapor 

p  =  constant  pressure 

8v  =  saturated  vapor 

t  =  at  temperature  T 

V  =  vaporization 

o  =  base  condition 
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PART  V 

GENERALIZED  THERMOOyNAMIG  AND  PHYSICAL 
PRO  PERTIES  OF  PETROLE  UM  TYPE  FUELS 


INTRODUCTION 

There  have  been  numerous  articles  in  the  literature  describing  various 
physiochemical  proper  ties  of  hydrocarbon  fuels.  For  this  information 
to  be  of  practical  use  to  industry  however,  it  was  desirable  to  have 
the  material  condensed  and  compiled  in  a  form  similar  to  a  handbook. 

The  data  presented  in  this  •  Part  V  fulfills  such  a  function. 

Much  of  the  material  contained  in  this  report  was  previously  issued  as 
Section  II  of  Engineering  Report  No.  206  prepared  by  the  West  Coast 
Laboratory  of  Thompson  Products,  Inc.,  Inglewood,  Calif.  ING.  ER.  206 
was  one  of  a  series  of  Interim  Reports  issued  under  United  States  Air 
Force  Contract  No.  AF  33(616)3729.  This  revised  Part  V# 
contains  empirically  derived  correlations  necessary  to  determine 
generalized  thermodynamic  and  physical  properties  of  petroleum  type 
fuels.  The  material  contained  in  Part  V  completes  the  study  of 

physiochemical  properties  of  selected  military  fuels. 

OBJECTIVE 

This  generalized  study  of  thermodynamic  and  physical  properties  of 
petroleum  type  fuels  was  made  for  one  basic  purpose: 

1)  The  determination  of  the  availability  of  generalized 
thermodynamic  information  on  hydrocarbon  fuels,  by 
means  of  a  literature  survey.  Any  reliable  information 
obtained  was  to  be  presented  as  part  of  this  report. 

SUMMARY 

All  reliable  information  that  could  be  obtained  related  to  generalized 
thermodynamic  and  physical  properties  of  petroleum  type  fuels  is 
presented  in  this  report.  A  compilation  of  charts  and  nomograms  has 
been  included.  These  charts  and  nomograms  were  obtained  from  an 
extensive  literature  survey.  They  present  what  is  believed  to  be  the 
most  accurate  correlations  of  the  physical  properties  of  hydrocarbon 
fuels  available.  In  some  cases,  the  original  charts  have  been  modified 
slightly  to  conform  to  the  limitations  of  the  fuel  types  being  investigated. 

I6l  . 
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DISCUSSION  OF  PHYSICAL  DATA 


The  charts  presented  in  this  report  are  a  collection  of  easily  used 
correlations,  covering  the  physical  properties  of  petroleum-type 
aircraft  fuels.  The  data*  collected  from  an  extensive  literature  survey, 
were  compared  and  analyzed.  A  comparison  of  this  data  with  available 
test  data  then  was  made.  In  cases  where  conflicting  information  was 
found,  the  source  which  best  approximates  the  average  has  been  pre¬ 
sented. 

In  some  instances  a  number  of  methods  for  predicting  properties  with 
varying  degrees  of  accuracy  are  presented.  In  each  case,  greater 
fundamental  knowledge  of  the  fuel  increases  ^he  accuracy  of  prediction. 

All  boiling  points  are  referenced  to  an  ASTM  distillation. 

The  characterization  factor  used  is  defined  as  »he  product  of  the  cube 
root  of  the  molal  average  boiling  temperature  in  degrees  Rankine  and 
the  reciprocal  of  the  specific  gravity  at  60°F. 

3 

K 


The  volumetric  average  boiling  point,  unless  otherwise  defined,  shall 
be  taken  as  the  arithmetic  mean  of  the  10,  30,  50,  70,  and  90  per  cent 
temperatures  of  an  ASTM  distillation. 

T  T  T  T  T 

VABP  =  10  30  h  50  70  +  90 

1 
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GENERAL  PHYSICAL  INFORMATION  AND  PROPERTIES 

DENSITY  -  GRAVITY  RELATIONS 
See  Figure  91 

Using  this  graph  the  specific  gravity  at  60*^  can  be  converted  to 
pounds  per  gallon,  pounds  per  cubic  foot,  or  degrees  APL 

VOLUMETRIC  EXPANSION  COEFFICIENT 
See  Figure  92 

Shown  in  this  chart  is  the  correction  factor  for  the  effect  of  temperature 
on  the  liquid  expansion  of  petroleum  fractions.  Jn  general,  the  useful 
range  of  this  chart  is  between  O’T'  and  300°F.  The  chart  is  also  useful 
in  determining  volume  changes  due  to  temperature  as  an  alternate  to 
the  use  of  density  relationships.  An  example  illustrating  the  use  of 
this  correction  factor  follows: 

Problem: 

Find  the  change  in  volume  due  to  a  temperature  change  of  75“F 
for  100  ft ^  of  fuel  which  has  a  40“  API  gravity. 

Answer: 

_4  3  3  ^ 

From  chart  oC  =  4.74x10  ft  /ft  ^ 

Change  in  Volume  A  V  =  x  A  T  x  V^ 

-4  3 

AV  =  4.74  x  10  X  75  x  100  =  3.  555  ft 

VARIATION  OF  VAPOR  PRESSURE  WITH  FUEL  TEMPERATURE 
See  Figure  93  '  j 

Vapor  pressure  as  a  function  of  volatility  factor  is  presented  on  this 
chart.  Volatility  factor  is  defined  as  the  sum  of  the  ASTM  volumetric 
average  boiling  pointland  the  ASTM  initial  boiling  point. 

MOLECULAR  WEIGHT  AS  A  FUNCTION  OF  BOILING  POINT  AND 

GR.WITY 

See  Figure  94 

The  plot  shown  in  this  figure  presents  a  moans  of  estimating  the 
average  molecular  weight  of  a  petroleum  fr.iction.  The  estimated 
deviation  is  plus  or  minus  ten  per  cent. 
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GENERALIZED  P-V-T  RELATIONS 
See  Figure  95 

The  usefulness  of  this  plot  cannot  be  overemphasized.  Nearly  ali  d<-vi- 
ations  from  ideality  can  be  computed  using  the  compressibility  factor 
indicated.  Pseudo-critical  temperature  and  pressure  should  be  used 
to  determine  reduced  temperature  and  pressure.  Informatuui  on  pscudo- 
critical  temperature  and  pressure  will  be  found  in  paragraphs  of  this 
Supplement  entitled  'Critical  Properties’.  Data  on  reduced  conditions 
greater  than  those  shown  in  Figure  95  can  be  obtained  from  the 
original  source.  (See  Reference  L  ) 

LIQUID  ENTHALPY  -  SPECIFX  HEAT  NOMOGRAM 
See  Figure  96 

If  two  conditions  are  known(temperature  and  specific  gravity  at  60“F), 
use  of  this  information  with  the  data  plotted  in  Figurc96 '^iH  permit 
determination  of  the  enthalpy  (above  0°F),  the  specific  heat,  and  the 
gravity  of  liquid  fractions.  The  maximum  error  estimated  for  this 
plot  is  plus  or  minus  five  per  cent. 

PROPERTIES  OF  PETROLEUM  FRACTIONS  NOMOGRAM 
See  Figure  97 

This  nomogram  is  quite  useful  as  a  rapid  means  of  determining  the 
following  properties  of  hydrocarbon  fractions: 

a.  Gravity 

b.  Carbon/  Hydrogen  weight  ratio 

c.  UOP  'K'  Characterization  factor 

d.  Aniline  point 

e.  Latent  heat  of  vaporization 

f.  Molecular  weight 

g.  Mean  average  boiling  point 

RELATIONSHIP  BETWEEN  VARIOUS  AVERAGE  BOILING  POINTS 
See  Figure  98 

A  knowledge  of  the  volumetric  average  boiling  temperature  is  required 
to  use  the  plot  presented  in  this  figure.  The  plot  shows  the  relationship 
of  the  other  average  temperatures  to  the  volumetric  average  boiling 
point. 
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HYDROGEN  CONTENT  OF  HYDROCARBONS 
See  Figure  ^ 

This  plot  presents  :  1)  a  means  of  estimating  the  percentage  of  hydrogen 
in  a  hydrocarbon  fraction  to  be  used  in  determining  heating  properties 
and,  Z)  the  solution  of  chemical  reaction  problems, 

APPROXIMATE  CHANGE  OF  SPECIFIC  GRAVITY  WITH  TEMPERATURE 
See  Figure  JOO 

This  plot  provides  information  concerning  the  variation  of  specific 
gravity  with  temperature  in  greater  detail  than  the  data  plotted  in 
Figure 


CRITICAL  PROPERTIES 


Numerous  correlations  have  been  based  on  the  theory  of  corresponding 
states;  i.  e.  ,  at  the  critical  conditions,  or  proportion  thereof,  most 
substances  behave  in  the  same  manner.  For  example,  the  physical 
properties  of  one  substance  at  conditions  having  half  the  value  of  the 
criticals  are  the  same  as  those  of  another  substance  at  the  corres¬ 
ponding  condition  (  relative  to  its  critical  condition). 

Reduced  Pressure  -  Ratio  of  actual  pressure  to  critical 

pressure  (both  absolute  values) 

Reduced  Temperature  -  Ratio  of  actual  temperature  to  critical 
__  temperature  (both  absolute  values) 

Mathematically  p  _  P 

r  *  P 

C 


T  * 
r 


T 

T 


In  the  case  of  mixtures  such  as  petroleum  fractions,  use  of  the  actual 
critical  properties  shows  poor  correlation  with  the  theory  of  corre¬ 
sponding  states.  To  improve  the  correlation,  a  pseudo-critical  state 
was  conceived  w^hich  can  be  related  to  the  true  critical  conditions. 
Thus,  in  the  case  of  petroleum  fractions,  the  reduced  conditions  are 
given  mathematically  as  follows; 
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r  P  P  =  pseudo  critical  pressure 

pc  pc 


T  s  pseudo  critical  temperature 
pc 

In  this  Section  plots  of  focal  temperatures  and  pressures  are  shown. 
Although  these  functions  have  no  theoretical  basis,  they  have  been 
derived  as  an  empirical  means  of  determining  the  equilibrium  vapor¬ 
ization  characteristics  of  petroleum  fractions.  The  focal  temperature 
and  pressure  represent  the  point  on  a  Cox  chart  to  which  the  various 
vaporization  percentage  lines  converge.  Knowledge  of  the  equilibrium 
vaporization  characteristics  at  some  pressure  plus  the  focal  point, 
allows  the  estimation  of  the  equilibrium  vaporization  characteristics 
at  other  pressure  levels. 

CRITICAL  PRESSURE 
See  Figure  101 

In  order  to  determine  the  critical  pressure,  a  knowledge  of  the  ten-to- 
ninety  per  cent  slope  of  the  ASTM  distillation  curve  along  with  the  API 
gravity  of  the  fraction  is  required.  In  the  example  show'n  on  the  plot, 
a  fuel  having  a  volume  average  boiling  point  of  276°F,  a  gravity  of  60® 
and  a  ten-to-ninety  per  cent  slope  of  3.  0  would  have  a  critical  pressure 
of  496  psia. 

PSEUDO-CRITICAL  PRESSURE 
See  Figure  102 

This  plot  is  straightforward,  requiring  a  knowledge  of  only  the  mean 
boiling  temperature  and  the  API  gravity. 


FOCAL  PRESSURE 
See  Figure  103 

Knowledge  of  the  critical  pressure  is  required  for  the  determination  of 
focal  pressure  from  this  plot.  The  value  obtained  from  the  ordinate 
is  the  difference  between  the  focal  and  critical  pressures,  thus 
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CRITICAL  TEMPERATURE 

See  Figure  104 


The  difference  between  the  critical  temperature  and  the  volumetric 
average  boiling  point  is  plotted  as  a  function  of  gravity  and  volumetric 
average  boiling  point. 

PSEUDO-CRITICAL  TEMPERATURE 

See  Figure  105 

This  plot  is  straightforward,  requiring  a  knowledge  of  only  the  molal 
average  temperature  and  the  API  gravity. 

FOCAL  TEMPERATURE 
See  Figure  106 

Knowledge  of  the  critical  temperature  is  required  for  the  determination 
of  focal  temperature  from  this  plot.  The  value  obtained  from  the  ordi¬ 
nate  is  the  difference  between  the  focal  temperature  and  the  critical 
temperature,  thus 


L  lOUID  SPECIFIC  HEAT 


The  values  of  specific  heat  given  in.  this  Part  V  are  shown  in 

three  plots.  The  difference  in  the  plots  is  the  degree  of  estimated 
accuracy.  If  K  factor  and  API  gravity  are  available,  the  plot  shown  in 
Figure  107  is  recommended. 

LIQUID  SPECIFIC  HEAT  AS  A  FUNCTION  OF  K  FACTOR  AND 
GRAVITY  AT  eo'F. 

See  Figure  107 

The  estimated  accuracy  of  this  plot  is  i  two-per-cent  for  petroleum 
fractions  similar  to  aviation  fuels. 
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APPROXIMATE  LIQUID  SPECIFIC  HEAT  AS  A  FUNCTION  OF  GRAVITX 
See  Figure  208 

This  plot  is  based  on  a  'rule  of  thumb'  equation 


and  is  a  ±  ten-per-cent  estimate  of  specific  heat  for  use  in  eases 
where  only  the  density  of  a  fraction  at  any  temperature  is  known. 


LIQUID  SPECIFIC  HEAT  AS  A  FUNCTION  OF  GRAVITY  AT  60“F 
See  Figure  209 

The  accuracy  of  this  plot  is  estimated  at  ±  five-per-cent  for  a  wide 
variety  of  petroleum  fractions.  The  plot  is  based  on  Cragoe's  equation^ 
Reference  2. 


VAPOR  SPECIFIC  HEAT 


There  appears  to  be  some  differences  between  various  investigators 
regarding  the  values  of  vapor  specific  heat.  The  specific  heat  plots 
presented  herein  represent  the  best  average  of  eight  correlations. 

VAPOR  SPECIFIC  HEAT  AT  ONE- ATMOSPHERE  PRESSURE 
AS  A  FUNCTION  OF  K  FACTOR  AND  GRAVITY  AT  60°F 
See  Figure  HO 

If  the  gravity  and  K  factor  information  is  available,  this  plot  should 
provide  data  with  an  accuracy  estimate  of  plus  or  minus  two  per  cent. 

VAPOR  SPECIFIC  HEAT  AT  ONE-ATMOSPHERE  PRESSURE 
AS  A  FUNCTION  OF  GRAVITY  AT  60‘F 
See  Figure  ill 

III  cases  where  only  gravity  is  known,  this  plot  can  be  used  with  a 
resulting  accuracy  approximating  plus  or  minus  five  per  cent. 
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PRESSURE  CORRECTION  FOR  MOLAR  JIEAT  CAPACITY  OP  GASES 
See  Figure  112 

The  effect  of  prcsaitre  on  the  specific  heaf  at  constant  pressure  is 
given  on  this  chart.  The  correction  factor  is  the  difference  between 
the  specific  heat  at  a  pressure  sufficiently  low  for  the  gas  to  act 
ideally  and  the  specific  heat  at  the  pressure  in  question  (with  both 
quantities  in  molar  units).  The  correction  is  added  to  values  obtained 

from  Figures  llbor  lllafter  making  the  proper  unit  conversion, 

GENERALIZED  HEAT  CAPACITY  DIFFERENCES 
See  Figure  213 

The  difference  between  the  specific  heat  at  constant  pressure  and 
constant  volume  is  shown  on  this  chart.  Here  again  the  quantities 
are  in  molar  units,  • 


HEAT  OF  VAPORIZATION 


The  term  heat  of  vaporization  is  quite  ambiguous  when  applied  to 
mixtures.  The  usual  definition  of  heat  of  vaporization  is  the  heat 
required  to  vaporize  a  unit  quantity  of  liquid  at  the  boiling  point.  This 
definition  requires  some  modification  for  use  with  a  mixture  having  a 
range  of  boiling  points  rather  than  a  single  boiling  point.’  The  values 
given  in  the  plots  in  this  section  are  the  heats  of  vaporization  assuming 
that  the  whole  material  boiled  at  the  mean-average-boiling-point  of  an 
ASTM  distillation. 

In  effect,  the  heat  of  vaporization,  as  referenced  in  this  discussion, 
is  the  heat  required  to  vaporize  a  unit  quantity  at  the  mean  average 
temperature  with  varying  pressure.  The  H  obtained  in  this  manner 
is  less  than  the  heat  required  to  vaporize  a  unit  quantity  at  one- 
atmosphere.  The  difference  between  these  H^  values  increases  as  the 
boiling  range  becomes  wider.  . 

Heat  of  vaporization  of  a  fuel  as  presented  here  can  be  viewed  also  as 
the  difference  in  enthalpy  at  constant  temperatu're  between  the  zero 
and  the  one-hundred  per  cent  vaporized  curves.  If  the  enthalpy  envel¬ 
ope  is  available,  the  heat  of  vaporization  at  constant  pressure  can  be 
determined  by  obtaining  the  difference  between  the  enthalpy  at  the  one- 
hundred-per-cent-vaporized  line  (at  the  end  point  temperature)  and  fl'c 
zero-per-cent-vaporized  line  (at  the  initial  boiling  point). 
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HEAT  OF  VAPORIZATION  AT  ONE- ATMOSPHERE  PRESSURE 
AND  MEAN  AVERAGE  BOILING  POINT 
See  Figure  114 

The  plot  can  be  used  if  the  mean  average  boiling  point  is  known  and  the 
molecular  weight  or  the  API  gravity.  It  is  recommended  that  the 
correlation  be  based  on  experimental  molecular  weights.  See  the  plot 
illustrated  in  Figure  97  for  an  additional  correlation. 

LATENT  HEAT  NEAR  THE  CRITICAL  TEMPERATURE 
See  Figure  I15 

At  the  critical  temperature,  the  latent  heat  of  vaporization  is  zero. 

The  latent  heat  of  vaporization  changes  abruptly  in  a  few  degrees  below 
the  critical  (approximately  50  BTU/lb.  in  50'^),  The  plot  presented 
covers  only  the  area  near  the  critical  temperature. 

In  order  to  use  this  plot,  a  knowledge  of  the  critical  temperature  and 
characterization  factor  is  required.  The  family  of  curves  shown  on  the 
plot  are  for  critical  temperatures  of  450,  500,  550,  600,  650,  700,  750, 
800,  850,  900,  950,  and  lOOO'F.  Interpolation  for  curves  having  critical 
temperatures  at  the  values  between  the  limits  of  400°F  t<  lOOO'F  can 
be  made  because  the  latent  heat  is  zero  at  the  critical  ti  mperature.  For 
example: 

The  heat  of  vaporization  at  650'F,  of  a  fuel  having 
a  K  factor  of  11.  0  and  a  critical  temperature  of 
675'F,  would  be  54.  5  BTU/lb. 

HEAT  OF  VAPORIZATION  CORRECTION  DUE  TO  TEMPERATURE 
See  Figure  II6 

This  correction  is  to  be  applied  to  the  heat  of  vaporization  data  obtained 
from  the  Nomogram  shown  in  Figure  97  or  the  Nomogram  shown 
in  Figure  114 

For  example: 

A  fuel  having  a  mean  average  temperature  of  400°F,  a 
critical  temperature  of  720“  and  a  molecular  weight  of  150, 
would  have  a  heat  of  vaporization  of  122.  5  BTU/lb.  (See 
the  Nomogram  shown  in  Figure  97  V  If  the  pressure  were 
increased  such  that,  at  the  new  pressure,  the  mean  average 
temperature  was  500“F,  then  the  latent  heat  would  have  to 
be  corrected.  From  the  plot  shown  in  Figure  7^6 
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T  -  T.  =  720  -  400  a  320“F.  T  -  T  a  720“  -  500"  =  Z20"F. 

C  k  ;  ■  C  . .  ■ 

The  heat  of  vaporization  at  500TF  is 

H  =  A  /  A  It  ^  @  1  atm 

H  =  (O.  86)  X  {122.  5)  s  105.  35  BTU/lb. 


ENTHALPY 


The  heat- content  property  of  petroleum  fractions  is  obtained  by  inte¬ 
grating  the  specific  heat  equation  over  the  temperature  range  desired. 
Normally  the  datum  is  taken  as  H  =  0  at  0"F. 

HEAT  CONTENT  OF  PETROLEUM  FRACTIONS 
See  Figure  117 

In  the  plot  shown  in  Figure  I17  the  enthalpy  of  the  liquid  is 
assumed  to  be  zero  at  0"F.  A  typical  example  is  illustrated  in  this 
plot  which  shows  the  method  for  determining  corrections  due  to  K 
factors  other  than  12.  0  and  pressures  above  one-atmosphere. 

HEAT  OF  COMBUSTION 


Both  gross  and  net  heats  of  combustion  plots  are  illustrated  in  the 
figures  applicable  to  this  topic.  The  difference  between  the  gross 
and  net  heat  is  the  heat  of  vaporization  of  the  water  formed  due  to 
combustion  of  the  hydrogen  present.  Only  the  net  heat  is  of  practical 
significance  in  normal  usage. 

GROSS  HEATS  OF  COMBUSTION 
See  Figure  II8 

Knowledge  of  the  API  gravity  and  K  factor  gives  the  gross  heating 
value. 

NET  HEATS  OF  COMBUSTION 
See  Figure  119 

The  plot  of  this  figure  is  similar  to  the  plot  shown  in  Figvirc  118 
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VISCOSITY 


Both  vapor  and  liquid  viacoaity  correlationa  are  included  in  the  dia- 
cuaaion  of  thia  topic.  Liquid  viscoaitiea  are  given  at  100*^  and  210*^, 
the  atandard  temperature  a  in  general  uae  by  the  petroleum  industry. 
For  values  at  other  temperatures,  graphical  extrapolation  or  interpol¬ 
ation  on  the  special  viscosity  chart  paper  of  the  Bureau  of  Standards 
is  recommended. 

ABSOLUTE  VISCOSITY  OF  HYDROCARBON  VAPORS 
AT  ONE-ATMOSPHERE 
See  Figure  120 

Vapor  viscosities  in  centipoise  are  shown  as  a  function  of  molecular 
weight  in  this  figure. 

KINEMATIC  VISCOSITY  OF  HYDROCARBON  VAPORS 
AT  ONE-ATMOSPHERE 
See  Figure  121 

In  this  figure,  vapor  viscosities  in  centisfokes  are  chosen  as  a 
function  of  molecular  weight. 

VISCOSITY  OF  GASES  AT  HIGH  PRESSURE 
See  figure  122 

The  effect  of  pressure  on  the  visco-  „ty  ratio  (i.  e. ,  ratio  of  viscosity 
at  pressures  above  one-atmosphere  to  that  at  one-atmosphere)  is 
shown  as  a  function  of  reduced  pressure  in  this  Figure.  In  correlation 
of  this  type,  the  pseudo-critical  pressure  should  be  used  for  the 
determination  of  reduced  pressure. 

RELATION  OF  VISCOSITY  AT  lOO'F  TO  CHARACTERIZATION 

FACTOR 

See  Figure  123 

The  average  boiling  point  shown  by  the  plot  in  this  Figure  is  the  molal 
average.  Use  of  this  plot  is  straightforward.  Knowledge  is  required 
of  the  API  gravity  and  either  the  molal  average  boiling  point  or  the  K 
factor;  viscosity  can  then  be  determined. 
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RELATION  OF  VISCOSITY  AT  210*^’  TO  CHARACTERIZATION 

FACTOR 

Sie  Figure  124 

The  average  boiling  point,  a«  illustrated  in  this  Figure,  is  the  molal 
average.  The  factors  applicable  to  the  plot  shown  in  Figure  123 
also  apply  to  this  plot.  Knowledge  is  required  of  the  API  gravity  and 
either  the  molal  average  boiling  point  or  the  K  factor;  viscosity  can 
then  be  determined. 
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APPENDIX  I 


The  correlations  found  in  this  Appendix  were  complied  and  cross 
checked  in  part  by  the  errpirical  relationships  outlined  in  the  following 
paragraphs, 

LIQUID  SPECIFIC  HEAT 

1.  Watson,  Nelson  and  Murphy.  Industrial  and  Engineering  Chemis¬ 
try,  Vol.  25.  1933,  p  880. 

C  =  0.  6811  -  0.  308/^+  T  (0.  000815  -  0.  000306/?;  (O.  055K  +  0.  35) 

P 

2.  Cragoe.  Bureau  of  Standards  Misc.  Publication  No.  97,  1929. 


:  =  - -  (0.  388  +  0.  00045T) 

P  0.  5  ' 


3.  Fortsch  and  Whitman.  Industrial  and  Engitieering  Chemistry. 
Vol.  18.  1926.  p  795. 


(210  -  (T  +  670) 
2030 
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7.  Goldstern.  Petroleum  Zeitschrift,  Vol.  32  No.  36,  1936.  pi. 

0.  402  +  0.  00045  (T  -  32) 

p  “  i 

8.  Mullins,  Journal  Institute  of  Petroleum.  Vol.  33.  1947,  p  44. 

Cp  »  0.  0009T  -  0.  499/^+  0,856 
where  T  is  in  Centigrade 

9.  Gaucher,  Industrial  and  Engineering  Chemistry.  Vol,  27. 

1935.  p  57. 

^  .  0.  623  -  0.  0187  (T)^  ^  0.  000355. (Jt  ^  o.  00045  (T  .  49i) 

P  .a  3/2  3 


where 

T  is  molal  average  in 

10.  Gaucher.  Industrial  and  Engineering  Chemistry.  Vol.  27. 
1935,  p  57. 


Cp  »  1, '323  -  0.2005K+  0.  0107K^  +  0.00045  (T  -  32) 

SPECITIC  HEAT  OF  VAPORS 

1.  Bahlke  and  Kay.  Industrial  and  Engineering  Chemistry.  Vol.  21. 
I929t  P.  942. 

I  4.0-^)  (T  t  570) 

S  6450 
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2.  Bahlke  and  Kay,  Induitrl*!  *nd  Enalnesrlng  Chemiitry,  VM.  21 

1929.  p  942. 

Cp  ■  (  4.  0 -  /^  )  X  (  0. 109  +  0.  00014T) 

where  T  la  In  “Centigrade. 

3.  Watson  and  Fallon.  National  Petroleum  News.  1944,  p  R372. 
Cp  =  (0.0450K  -  0.233)+ (0.440  +  0.  0177K  X  10"^T) - 

-6  2 
0.  153  X  10  T 

4.  Cragoe.  Bureau  of  Standards  Misc.  Publication  No.  97.  1929. 


0.  388  +  0.  00045  T  0.  09 


5.  Swelgert  and  Beardsley.  Experimental  Station  Bulletin  2, 
Georgia  School  Technological  Engineering.  1938. 


The  kerosene  vapors  having  molecular  weights  near  170  in  the 
temperature  range  >60  to  650°F. 


6.  Watson  and  Nelson.  Industrial  and  Engineering  Chemistry. 
Vol.  25.  1933.  p  880. 
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HEAT  OF  VAPORIZATION 

1<  KiitUkowtky.  Zeitung  Phyflik  and  Chemiiche.  Vol«  107«  1923. 
p  65, 

7. 58  T.  +  4.56  T  log  T 

H  ■  — '  — . . . 

MW 

where 

T^  •  mean  average  boiling  point. 

2,  Cragoe,  Bureau  of  Standarda  Miec.  Publication  No.  97.  1929. 

Hy  -  ~~~  (110.  9  -  0.09T)  X  — |— 

3.  Trouton  equation 

20.  5  T 

H  “  - - 

V  MW 

where 

T  is  th“  molal  av«>ragfl  boiling  point  in  Kelvin 


4.  Clausium  Clapeyron  equation 

Hy  ■  -If-  *  (V^  -  V  )  .  T 

where 


dP 

dT 


is  the  change  in  vapor  pressure  with  temperature; 


(  V  “  )  is  the  vapor  -  liquid  volume  difference  per  mole 

T  is  the  temperature  in  absolute  units. 
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5,  Goldatern,  Petroleum  Zeitschrift.  Vo),  32,  No.  36.  1936,  |>j. 

60  -  0.  05  (  T-32) 

«v  =  -  ^ 

HEAT  OF  COMBUSTION 

1.  Cogliano  and  Jessup.  Bureau  of  Standards  Report  2348.  March  It^ 
1953. 

6  2 

O  net  S  17,608  +  {0.2054  A  x/^)  -  7.  245  X  lo“  (A  x/^)  •  140(%s) 


where 


A  =  aniline  point  “F 
%8  =  per  cent  sulphur 

Cogliano  and  Jessup.  Bureau  of  Standards  Report  2348.  March 
18,  1953. 


O  net  =  22,  130  + 


2,  560 


/P  -  l.Si 


3.  ROJhberg  and  Jessup.  Industrial  and  Engineer! ng  Chemistry. 

Vol.  43.  1951.  p  981. 

Q  =  17,  944.  9  +  0.  1043  X  {.A  x/^  ) 

4.  Fein,  Nelson  and  She rman.  Industrial  and  Engineering  Chemistry. 

Vol.  45.  1953.  p  610. 

O  net  =  17,  145.9  +  0.  5155  (A  x -  0.  691 1  3  x  lo”*^  (A  x/5)'^T 
0.47772  X  10‘®  (A  )^  -  0.  1235  X  lo"^^  (A  x 

5.  Fein,  Nelson  and  Sherman.  Industrial  and  Engineering  Chemistry. 

Vol.  43.  1951.  pgai 


O  net  =  8,  505.  4  +  846,  aiK  +  U4,  92 4  0.  121865  -  9.  951K/® 
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CRITICAL  PROPERTIES 

1.  Roefls.  Journal  of  the  Institute  of  Petroleum  Technology.  Vol.  22 
1936,  p  665. 

T  critical  =  202.  7  +  1.  591a  -  6.  29  x  10““* 
where  a  =  /’x  (T  volumetric  average -t-  100) 

SYMBOLS 

C  =  Speci  fic  Heat  BTU/lb  'T*. 

H  =  Enthalpy  -  BTU/lb 

J  =  Mechanical  equivalent  of  heat 

K  =  Characterization  factor 

MW  =  Molecular  Weight 

P  =  Pressure  -  psia 

Q  =  Heat  of  Combustion  BTU/lb 

R  =  Gas  constant  -  10.73  psia,  ft^/'R,  lb-mole 

S  =  Entropy  -  BTU/lb  °F. 

T  =  Temperature  -  “F  or  “R 

V  =  Volume  -  ft^ 

Z  =  Compressibility  factor 

OC  =  Volume  residual 
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A  a  Increment 

X  »  Heat  of  vaporization  -  BTU/lb* 

P  ^  Density  at  60‘*F. 

SUBSCRIPTS  AND  SUPERSCRIPTS 


c  »  Critical 

k  a  Normal  mean  average  temperature 

L  =  Saturated  liquid 

m  =  Mean 

p  =  Constant  pressure 

r  a  Reduced  conditions 

8v  a  Saturated  vapor 

t  a  Total 

v  a  Vaporization 

0  to  20,  20  to  40,  etc.  =  0-20%,  20-40%  Fractions 

1  =  Initial  condition 

2  a  Final  condition 


Also  used  to  indicate  heat  of  vaporization. 


liiAEO  TR  60-767 


180 


REFERENCES 


The  source  of  plots  Illustrated  in  part  V  Figures  are  listed 
below. 

1.  Maxwell,  J,  B.  Data  Book  on  Hydrocarbons,  First  Edition, 
D.  Van  Noatrand  Company,  Inc.  ,  Princeton,  New  Jersey, 
March,  1951.  pp  21,  41,  72,  73,  148,  174,  175,  177. 

2.  Cragoe,  C.  S.  Bureau  of  Standards,  Misc.  Publication 
No.  97,  Nov.  9,  1929. 

3.  Barnett,  Henry  C.  and  Hibbard,  Robert  R,  Properties  of 
Aircraft  Fuels,  National  Advisory  Committee  for  Aero¬ 
nautics  Technical  Note  3276,  Lewis  Flight  Propulsion 
Laboratory,  Cleveland,  Ohio.  August  1956.  p  78. 

4.  Bureau  of  Standards  Data, 

5.  Findl,  E.  and  Edwards,  H.  Vaporization  Cycle  and  Thermo¬ 
dynamic  Properties  of  Jp-4,  RJ-1  and  JP-6-H.  INC.  ER234. 
Study  of  Physiochemlcal  Properties  of  Selected  Military 
Fuels.  Thompson  Products,  Inc.  Part  V,  Supplement  V. 
December  19-58. 

6.  Winn,  F.  W.  Physical  Properties  by  Nomogram.  Petroleum 
Refiner.  Feb.  1957.  p  157. 

7.  Watson,  K.  M.  and  Nelson,  W.  L.  Industrial  and 
Engineering  Chemistry  Vol.  25.1933,  p880;  Vol.  29,  1937, 
p  1408. 

8.  Hougen,  O.  A.  and  Watson,  K.  M.  Chemical  Process  Prin- 
ciples  Charts.  John  Wiley  and  SonsT^Inc. ,  New  York,  N.  Y. 
April  1956. 

9.  Nelson,  W.  L.  _OiI  and  Gas  Journal.  January  27,  1938. 

P  184. 

10.  Edmister,  \V.  C.  and  Pollock,  D.  H.  Chemical  Enyineerin^ 
Progress.  Vol.  44  No.  12.  pp9U,  924,  925. 

1 1.  Holcomb  and  Brown,  Industrial  and  Engir.eering  Chemistry. 
Vol.  24,  1932,  pZlO;  Vol.  18,  J  iZb,  p  795. 

WADD  TR  60-767  181 


12.  Edmister,  W.  C.  The  Petroleum  Engineer.  Dec.  1950. 
p  C-16. 

13.  Edmifiter,  W.  C.  Petroleum  Refint-r.  Nr*v.  1948.  p  C-16, 

14.  Nelaon,  W,  L.  Oil  and  Gas  Journal.  July  26,  1947.  p  254. 

15.  Watson,  K.  M.  Industrial  and  Eagineering  Chemistry.  Vol. 
35.  1943,  p  398. 

16.  Fine,  R.  Sy  Wilson,  H.  I.  and  Sherman,  J.  Net  Heat  of 
Combustion  of  Petroleum  Hydrocarbons.  Industrial  and 
Engineering  Chemistry,  Vol.  45,  1953,  p  610. 


17.  Bell,  H.  S.  American  Petroleum  Refining.  Copyright  1945, 
Third  Edition,  D.  Van  Nostrand  Company,  Inc.  Princeton, 
New  Jersey,  pp  54,  55. 

18.  Rothenberg,  A.  J. ,  Edwards,  H. ,  Woodell,  E.  W.  ,  Volatility 
Characteristics  of  Eight  Types  of  Turbine  Fuel  at  Tempera¬ 
tures  to  450°F.  ,  West  Coast  Research  and  Development 
Laboratory,  Report  No.  INC.  ER.  176,  Thompson  Products, 
Inc.,  June  15,  1954. 


Figure  91  Plot  o£  Denelty-Crerlty  ReUtlone 


M  ^ert^M* 


50 


too  150  200  250  300  400  500  600  700  800  90C 

TEMPERATURE  ®F 


yirur«»  93  -  Variation  of  Vapor  PT<*3sur>*  vlth  Fu^l 


Moltculor  Wtight  Moltculor  WtighI 


r  Relation  of  Hydrocarbon  Vapors  (Mole'cular  Wt.  Equal  to  or  Greater  than  16) 


TEMPCRATURE  OF 


640 


A  STOAiti^rr  Line  Conhectjmq  any  two  Pbopertics 
GIVES  ALL  OTHER  BELATED  PROPERTIES 


600 


(R-f.  5 

Fli;ur»  96  -  Prcp»rti»3  of  liquid  P<>tro\f'un  FraoHlons  Norropra? 
VAfD  TR  60-767 


168 
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Figure  109  •  Liquid  Specific  Heats  of  Petroleum  fractions 
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